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ABSTRACT 


The replenishment at sea (RAS) maneuver is studied in 
detail fcr heading and speed control. Design OL 
purposefuliy ncnlinear control laws is accomplished fcr the 
Mariner hull using the linearized equations of amcticn in 
three degrees cf freedon. Extensive use of lew order 
modeling and optimal control theory was. made. Procedure 
steps are presented in detail to facilitate redesign for 
other ship tyres. The results are verified using DSL 
Simulaticn fcr a number of possible RAS scenarios. The 
ccntrol systems are also tested in a sea state to insure 


preper cperation in the presence of external perturbations. 
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I. INTRODUCTION 


The advent of the digital computer as standard 
equipment cn toard virtually all modern Naval ships has 
opened the field of Digital Computer Contrcel in almost all 
aspects cf ship life. The computer has been a viable asset 
in fire ccrtrol systems for years and has been used 
extensively for aids to ship maneuvering control in the form 
es NTDS (Naval Tactical Data System) readouts. The 
declining ccsts cf general and special purpose computers has 
made their Loception as a Manpower replacemert or 
augmentation a reality. Their high speed and accuracy can 
Make them fperform functions with much greater safety than 
previously attainakle with time proven fand sometimes time 


weary) “seaman's eye." 


This then is the basis for this thesis; a study of the 
types of maneuvers that can be handled more accurately and 


safely than fresentiy being accomplished. 


One area of study is the total Replenishment At Sea 
(RAS) protlem including both course and speed control for 
tke apprceach and alongside phases. This situaticn has 
always been cne of extreme danger due tc the collision 
Fotentials invclved. However, other dangers are jnvclved in 
the on deck evclutions when the ships are not kept at a 
fairly cecnstant distance. Sudden violent maneuvers may 
cause the stress on the lines between shifs to increase 
enough tc cause the lines to part. The reality of this 
danger is readily apparent to anyone who has ever seen a 
Marila line cr steel cable part or a kingpost shackle break 
or a kingfost suddenly rend under these extreme stresses. A 


system which will minimize these potential dangers is well 
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werth investigation. 


Of course with a digital computer, the algorithm for 
RAS can ke modified or replaced to enable its use as a 
Maneuvering contrcl device for other situations such as 
formation steaming or single ship navigation transit 


eentrol. 
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IIT. MODELING 


A. MARINER CYNAMICS 


In the ccnception of this thesis, realistic models of 
modern destrcyer hull configurations were sought. This 
search proved fruitless. The hydrodynamic coefficients for 
present cay destroyers are not currently available. 
However, scme naval and civilian research is presently being 


conducted to cktain these coefficients. 


A cecmplete set of these coefficients is necessary for 
any maneuvering control system design. A hull configuration 
which has teen under continual study with well defined and 
verified hydrodynamic coefficients was chosenf{!]j., This hull 


is commonly referred to as the "Mariner" hull. 


The development of the equations of motion in six and 
three degrees of freedom have been well documented{*]. The 
mcdel used fer this thesis is the eguations of moticn in 
three degrees of freedom linearized with second order and 
higher terms eliminated. These equaticns are characterized 
by dependency on small perturbations about a specific 
Operating point. The maneuvers experienced in the fcllowing 
charters dc not entirely meet this criterion. The 
inadeguacy and shcrtcomings of this model are of little 
ccnsegquence Lecause relevant hydrodynamic caefficients are 
not available, and the methods presented can be applied to 
any ship tyre. 


The development of the model is readily available to 


the interested reader [3]; only a summary of the equations 


Us 





and theix ccrresfonding hydrodynamic GEetricientes are 


presented kere. 
The equations of motion used are as rollows: 


X.-2)U+X (U-u ) +X d=0 

( u ms o? d 

(Yo-u) V+Y V#(Y -m)R+Y.R+Y d=0 
V v 1B es d 


(No-I )F4N R4#N.V+N V+N ~d=0 
1a) a4 e v Vv d 


The direction and sense of the terms in the akove 
equations are shown in figure II-1. Letting: 


a =E-Y. 
11 Vv 
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Figure II-l 


Direction and Sense of Hydrodynamic Terms 


Ze 





The equaticns cf motion can Le written as: 


a Atb Atc Ata Btb Btc B=IF1 
“eed ite 21 #21 «21 


a Atb A+tc Ata Btb Btc B=IF2 
1D (55s i ye ey: 


a Ctb Cte C=IF3 
33 33 


33 
ors 
| tt 
am Ata .-=i4 
Zz 
rT tt 
a Ata E=I2 
2 
t 
a _C=I3 
where; 
IF1=-Y ed 
IF2= eq 


22 











re3e- f soy ane 
u 


and: 


I1=-—k &-c A-b B-c BtIF1 
7 4) a & 


I2=-t ~B-c Aw-b B-c B+tIF2 
one 22 


By sclving this system of equations, the following 


relaticnshifs are established: 


A=(a@ a =Ti-a_a [{2)/DEL 


we 25 es 21 


B=(a a I2-a a 11)/DEL 
1133 Bea pz 


DEL=e4 (2) el eh 
See 11 22 12 21 


which yields the solution: 


: tt t 1 
B dt=9 [ cB (a) +f Bdt jdt 
t> 2 to ty 


The space cocrdinate system is defined as fcllows: 
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Y=UeSING+Vecosyp 


X=U «COSP-VeSIND 


where: 


€ 
X=X “f X dt 
0 te 
€ 
Y=Y f 
OS, 


Table II-7 summarizes the symbols and nemenclature used 


<j o 


dt 


in the fcregcing abbreviated solution of motion in three 
degrees or freedom. The applicable characteristics of the 
Mariner hull are presented in table II-2 with the 
corresponding nondimensionalized hydrodynamic coefficients 
and the DSL computer program variable names delineated in 
maebe IT-s. | 


Computer Program #1 is the basic DSL program that was 
developed frcem these equations of motion. This prograkz uses 
two ships to illustrate the turning characteristics cf the 
Mariner hull fcr various rudder commands. Figure [I-2 shows 
the difference between a step mcdel and a ramp nodel rudder 
in a gecgraphic flot. Figure II-3 shows the correspcnding 
difference in yaw. 
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TABLE II-1 


SYMBOLS AND NOMENCLATURE 


derivative 


respect to 


derivative 


respect to 


derivative 


transverse 


derivative 


transverse 


derivative 


tc angular 


derivative 


to angular 


derivative 


Ep ab eS SSS OS Se 


Of longitudinal force With 


longitudinal acceleration JU 


of longitudinal force 


longitudinal velocity 


of lateral force 


acceleration V 


of lateral force 
velocity V 


of lateral force 


acceleration R 


of lateral force 


velocity R 


cf lateral force 


to rudder angle d 


with 


with 


with 


with 


with 


with 
U 


respect to 


respect to 


respect 


respect 


respect 


derivative of yaw moment with respect to 


transverse 


acceleration V 


derivative of yaw moment with respect 


tc transverse velocity V 
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TABLE II-1 (cont.) 


SYMBOLS AND NOMENCLATURE 





Definition 


derivative of yaw moment with respect to 


angular acceleration R 


derivative of yaw moment with respect tc 
angular velocity R 


derivative of yaw moment with respect to 
rudder angle d 


yaw angle acceleration 
yaw angle velocity 


initial velocity of origin cf body axes 


relative to fluid 


transverse acceleration of ship axes 
relative to fluid 


transverse velocity of orign of ship axes 
relative to fluid 


hydrodynamic longitudinal fcrce 


hydrodynamic lateral force 


longitudinal acceleration cf Ship axes 
relative to fiuid 
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TABLE Li’ (COnt .} 


SYMBOLS AND NOMENCLATURE 


Symbcl Definition 
U 


ee velocity of ship axes relative 


yaw angle 
actual time 
nondimensionalized time 


longitudinal distance that the ship center 


of gravity is forward of the ships axes 


longitudinal velocity of ship axes relative 


tc fluid (cperating pcint) 


TABLE £1—2 


SHanaG Lent Slees OF MARINER-TYPE STUDY SHIP 


Length, ft. B21 
Beam, £t 1600 
Draft, ft 29.75 
Displacement, tons 16,500. 
Block coefficient, e. 0.26 
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DAB io thas 


NCNDIMENSICNAL HYDRODYNAMIC COEFFICIENTS 


Computer Program Nondimensional 
Coeficient Variable Name Value 
(X.-b) MXUD -0.0085 
u 
X XU -0.0012 
u 
ye yO -~0.01243 
v 
(Y.—t) MYVD SUS) 
Vv 
(Y -mu ) MYR =U5005 7 
(Y.-ex ) YRD —0. 00027 
g 
3 YDELR +0.0027 
N NV “0.00351 
v 
N. NVD -~0.000197 
v 
(N -mx u ) NR -0.00227 
(N.-I ) IZNRD -0.00068 
an Zz. 
Me NDELR =-0.00126 
X eer -~0.0000462 
n 
ss a ~0.0000052 
n 
N --- +0.0000026 
n 
X XDELR 0.0 
d 
NOTE: x =C€.0 
g 


28 





Tapers 1 CGOnt.) 


NCNCIMENSIONAL HYDRODYNAMIC COEFFICIENTS 


Values Eased on the following operating point: 
a (15 Kts) 


2g 
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Bae BRUDDER RESPONSE 


The previous section indicates a marked difference in 
behavicr retween step and ramp rudder models. This prempted 
an investigation into realistic rudder modeling which would 
fulfill tte reguirements of limit stops and maximum rcudder 
rate. ; 

NSRCOC{*]} has modeled the rudder of the DD-931 Class 
Destroyer. The basics of this model are presented in the 
block diagram of figure II-4. 


The first limiter models the rudder stcfs which fcr the 
Mariner are +30 degrees. The second limiter models the 
proportional band of a variable-displacement pump by 
ligiting its Maximum percent stroke. The limits for this 


nonlinear element have been found to be +7 degrees. 
The transfer function (K ys) acceptS an input error 
9g 


Signal of up to 7 degrees, converts it to a rudder rate, and 


integrates tke rudder rate to obtain rudder angle. Letting: 
qd = Maximum rudder rate (2.0 degrees/sec) 
D 


d = Maximum error input (7.0 degrees) 
emax 


The system gain can Le defined as: 


K =4a@y 
g - emax 


2507 120 
Q0.2€5714 /sec 


5)% 
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Tc convert this model to the required nondimensionalized 


form, the fcllcwing manipulation is required: 


K §*' = K eL/u 
g g 1 
= 0.285714e527.8/(15¢1.689) 
= 5.95224 
wheres; 


L = shif length 


a = cpéerating point speed (15 Ktse1.689 ft/sec/Kt) 


Computer Program #2 is the DSL program which rodels 
this system. The curves of figure II-5 and II-6 exhibit the 
responses cf various step rudder commands. These are 


tabulated and cress referenced in table II-4. 


These responses show the characteristics of a realistic 
rudder in that the rudder is never allowed tec slam intc the 
stcps. They exhibit the time delay between command and 
response which is a function of the rate of response (2.0 
degrees/sec). A centrol system design with this scheme is a 
much more difficult problem than one with an idealized 
rudder (step response) because the entire rudder ccntrol 


system beccmes quite nonlinear. 
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TABLE II-4 


RUDDER COMMAND AND RESPONSE 


Rudder Senet at 
Figure Curve Command (deg) Condition(deqg) 


r—5 1 +30.0 -30.0 
ri—5 fs +25.0 SPANO 
i= > 3 +20.0 PG, 
ii 5 4 +#15.0 = 4h) 610, 
it) 5 +10.0 =O 0 
i= 5 6 tao) = 5-0 
iG 1 +3C.0 O50 
i= 6 2 25.0 0.0 
i= 3 +20.0 0.0 
T= 6 4 tlle G 0.0 
TI-6 5 +10.0 0.0 
© 6 = 530 0.0 
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C. ENGINE RESEONSE 


Figure II-7 pertrays a complex model of a gas turbine 
Peogrulsicn filant{ *j[°]. This model cantains the elements 
required for a complete dynamic study of the systen. lo) 
the purpese ci this thesis, such a complicated model is not 
required if the overall input-output relationship can be 
estaklisted. 


Reference 5 establishes ane OUtPUT speed (0) 
relaticnship for a step input of desired speed (Ud) and is 
redrawn as figure II-8. The relationshifp aprears to be that 
of a first or seccnd order system with a time delay. 


The system equations for a first order approxization 
With a tine delay say be written as: 


SPLTIN = K*SPDDES*¥e 

SEDERR (SEDIN-SEDCUT) *G 

SPLOUT = f severe dt 
Which yeilds the transfer function: 


S 
SELGUT(s) G*K*e 


SPLIIN(s) $+G 


Which is bleck diagrammed in figure II-9. 





SPDDES a e _ SPDouT 


| Figure II-9 
Propulsion Plant Low Order Model Block Diagram 
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From figure II-8, the time delay, system gain, and time 
ccenstant can be estimated as: 
T = 4.88 sec 
0.9877 
g2 092 


QQ 
iT) 


Computer Program #3 was used to okttain the ster input 
response. The original complex system output and the low 
order approximaticn are compared in figure II-10. As 
indicated in this figure, the two responses are very close. 
Considering the linearized approximations made in the 
equations cf mction, this response 1s accurate encugh for 
system study use and is used as the model fcr speed ccntrol 
in chapter III. 


Similar methods may be used to oktain simplified 
low-order acdels for other. high-order propulsion systems now 
in use (e.g. pressure fired boiler systems, 1200 1b. 
systems, etc.). They may not, however, simplify to a first 
order afrreximaticn suitable for system study. A method of 
computer determination of low-order models of high-crder 
Systems is ccntained in ref. 7 and may be mechanized for 
this purpese. 
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D. EXTERNAL FORCES 


The modeling of ship dynamics cannct be complete 
without the introduction of external forces which perturb 
its respenses. These forces are caused by many factors and 
some are mcre relevant than others in the scope of this 
thesis. The two that are considered can cause substantial 
perturbations that must be modeled and eventually acccunted 


for in tke ccatrcl system design. 


Whenever two ships operate in close proximity (less 
.than 250 feet), suction and pressure forces between hulls 
are present. Studies have peen conducted on the Mariner 
hulif[!] which have produced data for ccnstruction of a 
fanily of curves fcr two ships passing on the same heading. 
No data has teen found for the cases of two ships not cn the 
same heading. Other restrictions on the work presented in 
ref. 1 are that the ships are of the same tyfe and of 
Similar hull ratios. 


Interactive effects vary as the square of speed. 
However, this 1s only true if both ships are at the same 
speeds. The interaction modification factor is based cn the 
normalizec speed of 15 kts. This factor can thus be written 
ass 


2 
SPDE = CDOT! 


Exact effects on the interacticn forces and moments in the 
Situation where the ship's speeds differ are not available. 
This is inconvenient since during the approach phase, the 


nciImalized speed of the approach ship (ship 5B) can fe as 
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high as 1.5. If the effect on ship B is as stated above but 
With its cwn speed causing the interaction modification, the 
interactive forces and moments can be 2.25 times greater 


than withcut speed considerations. 


Withcut the ability to pin down this relationship, it 
was ignored in the development cf the contrcl laws presented 
in this thesis. Appendix C was written with the expressed 
intent of illustrating the effect ene modifying the 
interactive forces and moments to the extremum case 
menticned akcve. It must be realized that this case is not 
ccnsidered likely in that it is felt that the interactive 
forces and mcméents modification cn ship B are more apt to be 
caused Ey the speed of ship A. If this iS so, since ship A 
is kept at a constant 15 kts., the interactions need ano 
mcdificaticn fcr speed consideration in this thesis. 


Reference 1 also gives a method of modifyince the 
interactive forces and moments based on different ship 
lengths. FCr ease of computation and graphical 
presentation, the two ships were considered of equal 
lengths. Tc modify this to ships of dissimilar lengths, the 
resulting hydrcdynamic derivatives must be mcdified as shown 


in ref. 1 (alsc shown in appendix Cc). 


Since neo closed form expression existed for these 
forces, tke fagily of curves reproduced as figures II-11 and 
II-12 were guantized in subroutine SLOPES (an adaptation of 
the subroutine of the same name from ref. Vi)” (Appendix Cc 
contains a curve fitted subroutine that was compiled after 
tke completicn of the research on this thesis. It was not 
used for any design qc simulation runs except for those 
presented in that same appendix.) An interpolation 
algorithm is used to approximate the intermediate values 
Fetween guantized values and between the curves of the 


fawaly. 
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Lateral Distance DY (Ft) 


-400 -200 0 200 400 


Note: To convert to normalized lateral distance DY - 
Divide by the ship length L. 
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Family of Interactive Y force Curves 
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* 
Note: To convert to normalized lateral distance DY - 
Divide by the ship length L. 
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Family of Interactive N Moment Curves 
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The main purpcse of this subroutine is to compute the 
interactive forces between ships in the replenishment at sea 
Situation anc output the values for perturbation of the 
centrol ship only. Even though both ships are affected by 
‘these perturkaticns, a one ship control system which is 
effective regardless Oe the other ships moticn is 
considered. Ccnseguently, the interactive forces on the 


seccnd ship are ignored. 


Subroutine SLCPES is contained in appendix A. Figure 
IiI-13 is a geographic plot of the two ships rfassing at 105.6 
feet with their rudders amidships (0 degrees). Figure II-14 
and figure II~15 show the magnitude of the lateral force (Y 
ferce) and rctational moment (N moment) of the reference 
Ship on the ship making its approach (control ship). The 
reference ship is at 15 kts. and the apprceach ship is at 
Pees KtsS. The control ship starts its approach 5 ship 
dengths (2639.0 feet) astern and 0.4 ship lengths (211.12 
feet) lateraily displaced. The most graphic portrait cf the 
effects cf these fcrces and moments appears in the yaw 
Changes which are presented in figure II-16. From these 
figures it beccmes readily apparent that these perturbations 
cause viclert motions of the ship which must be acccunted 
for in any ceontrcl system development. Throughout the 
development of stich a control system in chapter III, these 
forces and mwscments are considered inherent in the model for 
RAS ccntrcl. 
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2. Waves 


The nodeling of sea state in the form cf waves and wave 
interacticns has cccupied the time of many researchers 
(Sj[?]. The exact formulaticn of waves will act be 
accomplished in this thesis. Since the main concern here 
Will be toc test the control scheme developed in chapter III, 
a much sigzplified wave generator can be used. To introduce 
the required experimental perturbations cn the designed 
control system a periodic wave system with a fundamental 
frequency ard its second harmonic is used. Some small 
Candom wave properties are introduced that ride on these two 
sinusoids. A simple expression of this combined wave can be 


written as: 


2 
W = WFesin(WE)+(ELeWF /WL) esin(2¢WE) +WFeWRVesin (WE) 


") denctes the Wave 

WF ais the Wave Force 

WE is the Wave Encounter radian Frequency 
WL is the Wave Length 

WRV is the Wave Random Variable 

PI is =.1415926 


With this wave as a basis, a method cf acdeling this in 
the dynamic environment of the total ship simulation was 
defined. The modeling includes the intrcduction of this 
wave intc the three degree of freedom eguaticns of mction. 
To accomplish this a set of defining relationships were 
develcped. First the general wave directicn is input to 
establish the direction of the wave encounter on the ship. 
If the ship direction is YAWDP2 and the wave directicn is 


HD, the expected wave direction is defined as: 


EWC = WD-YAWDP2 
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Next the wave enccunter frequency (radian frequency) can be. 
established with knowledge of the ships nermalized true 
speed (CDCT2), wave length (WL) and normalized wave velocity 


(WV). The wave encounter frequency (WEF) is then: 
WEF = zePIe (CDOT2+WVecos (EWD))/WL 


The total wave encounter (WE) is ncthing more than the 
wave enccunter frequency (WEF) times time. This gives the 
wave enccunter radian frequency required in the simple 


expressicn fcr the combined wave previcusly shown. 


This dces not ccmplete the task, since the individual 
wave forces of each degree of freedom must te derived for 
this general wave expression, namely the ccmponents of WF. 
Again a much simplified version of the more complex real 
life wave forces were used. The X and Y forces are 
considered first. These can te modeled as cosine and Sine 


functicns cf the expected wave direction (EWC) such that: 


WE X 
WEY 


WF ecos (EWD) 
WPesin (EWD) 


where WF is the total wave force of the enccuntered wave. 


The rotational N forces are a little more difficult to 
establish. Ey ccnsidering that no rotational forces are 
created Ey a wave directly on the bow or stern, or directly 
cff the teanz, and that it is maximum when the wave is at 45 
degrees off the bow Or stern, a much Simplified 
apEroximaticn is developed. Realizing that this method is 


very crude, the N force can be written as: 


WEN = WFesin (2eEWD) 


To add mere creditability to the wave mcdel, a random 


ug 


variable is added to the wave force at the waves fundamental 
freguency. A gaussian (normal) distributicn was chosen with 
a zero mean and a standard deviation of one-tenth the 
Maximum allowarle force of WF. A zero mean signifies that 
the expected amplitude of the random wave is 0.0, while the 
standard deviaticn Signifies that 68% of the random waves 
will have agplitudes less than one-tenth of the maximun 
allcwable ferce of WF. Also, 94% will have amplitudes less 
than one-half the maximum allowakle force of WF. This small 
added perturkation allows for verification of the model 
Simulaticn wath a stochastic force, which in turn adds 


creditability to the developed ccntrol systems. 


What remains is to define the total wave force (WF). 
It is impertant not to fall into a common Simulation pitfall 
Which inevitably causes unneeded design changes. A sea 
state does nct increase at an infinite rate. It therefore 
is incorrect to start a simulation with initial conditions 
set for calm séa and immediately introduce a high sea state 
perturbaticn. The initial large perturbaticn transient can 
give results that are not only unrealistic, but can cause 
the model and control system to produ3e unstable results. 
This is especially true in this case since the linearized 
(stfall fperturkaticns about an operating pceint) equaticns of 


mcticn are uwsed. 


With this in mind, a ramp feed in of the wave force 
with a limiter at the desired maximum wave ferce (WFMA) was 
used. The slcpe of the ramp was established to impart 
Minimum initial transients, yet increase the wave force to 
an acceptable testing level within the time frames cf the 
Simulaticnos cf chapter III. The slope is designed such that 
the maximum wave force is reached in 94.815 seconds actual 


Simulaticn time (4.548 seconds problem time). 


Computer program #4 was used to verify the wave action 
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model. Table II-5 on page 56 indicates the figures frecduced 
and changes in input wave length (WL) and wave direction 
(WD) for each run. The input parameters that were constant 
for all runs are tabulated below: 

YARDE2 = 0.0 

CEGEZ = 1.5 


” 
WS = 5.0 
WFMA = 0.1137 


x 
NCTE; WS is the unnormalized wave speed. Conversicn to 
ncrgalized wave velccity is: 
WY = WSY15.0 


Intrcduction of the wave forces is accomplished by 
multiplying the established wave forces by the rudder 
hydrodynasic coefficients for the individual reference 
directions. This effectively scales the wave forces tc the 
ship model tkéing used. The wave rorce result is coded in 


tke ship simulation program as follows: 


IF12 = KA1eD2+YY2+KAIOWY 
IF2Z = KELTODZ+YN2+KBITOWN 
IF32 = KC1eD2+NC2+KC1 oWX 


Detailed results of the wave farce effects are given in 
Chapter III and will not be dealt with here. 
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* 
NOTE: WL is given in ship lengths 
WD iS given in degrees 
Curve numbers of all runs corresponding to 
wave force components are: 


Force Curve 
WX ag 
WY 2 
WN S 
Table II-5 


Wave Simulation Listing 
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PEt REPLENISHMENT AT SEA 


Aj HEADING CCNTROL 


fee Contrel Chcice 





Many studies involving replenishment at sea (RAS) have 
treated the preblem asa multivariable system{ 19 j[t1tj}[12] 
{43}. Acadetically, there is nething wreng with this 
approach. Heowever, aS a practical system it leaves much to 
be desired. The key drawback in the multivariable system is 
the inescapakle dependency cn a command and control link 
between tke replenishment ships. The unreliability of UHF 
communications at these close distances is a much 
experienced phenonenom. It is felt that any knowledgeable 
cCamanding cfficer would not entrust the safety of his ship 
to such a guestionable link. An alternate method which is 
described kere is a modern extension of the long trusted 
"seaman's eye" ccncept, Where the senscrs and ccntrol 
devices mest be self contained on the Ship making the 
@meecach (tereafter referred tc as the receiving ship or 
eae 8B). 


In all fresent day RAS operations, the ship on which 
the apprcach is being made (hereafter referred to as the 
supplying shir or ship A), must maintain the replerishment 
course and speed. The receiving ship accomplishes the 
Maneuvers tc Baintain station relative tc the supflying 


Simic 


The parameters which are presently measured "visually" 
are relative position (in both the X and the Y directions), 
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relative heac (usually in reference to crdered replenishment 
course), and relative motion in the X direction (fcr speed 
matching). ‘These parameters are usually visualized by the 
ccnning cfficer who in turn gives corrective orders to the 
helzsman. The helusman must then translate these verbal 
orders into rudder and speed ccmmands through the helm and 
lee helm ccnsecles. The accuracy of the execution of the 
conning cfficer‘s ocrders 1S extremely dependent cn the 
ability cf the helmsman and throttleman. This system can be 
guite effective, and it can also be quite disastrous. MThis 
fast reacting and constantly changing environment lends 
itself tc Ereakdcwn in communications and manifests the 
inability of scme individuals to cope with the required 


critical ran-sachine interfaces. 


To eliminate these problems, present state of the art 
digital ccmruters and sensors are available for immediate 
itplementaticn of a completely automatic ship Centro. 
system. Such a control system may be installed on 
individual ships and be used for RAS without the reguirement 
of having the matching installation on the cther ship cf the 


heckup (another drawback of the multivariable approach). 


One cf the many pitfalls that may be encountered in 
digital simulation is the reality of the parameters that are 
measurable ir the real world situation as opposed to those 
that are incidently available in the simulation. With this 
fact as a keynote, Subroutine RBMEAS (Range and #éearing 
MEASurement) was developed. This subroutine, as listed in 
appendix Aa, defines the forward (FWD) and after (AFT) 
relative and true bearings, and ranges from the receiving 
Ship to similar points on the supplying ship. Figure III-1 


delineates tke terms used in the subroutine. The SDFn terms 
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are used tc position the bow and stern sensors and 
reflectors cn ship B and A, respectively, in geographic 
cocrdinates as a function of ships head. The FWD distance 
on the xX cccrdinate is ADPX and the Y cocrdinate is ADFY. 
Sizilarly, tte AFT distances are ADAX and ADAY. R1 and R2 
are the FWE and AFT ranges measured by a highly accurate 
ranging device installed on ship 8B. This same fLCanging 
device, if properly provided with a pinpcint reflector on 
the supply -ship(ship A), will give accurate relative 
bearings FWD and AFT., B1 and B2 respectively. The distance 
between senscrs may be varied, but as a rule should te kept 
as far apart as pessible to allow maximum sensitivity. The 
distance used in this thesis is 1.0 (one ship length), and 
the distances were considered the same for beth ships. This 
is net a necéssary condition and may be changed to suit the 


Situaticn. 


Subreutine RBKEAS assumes highly accurate sensors in 
both rance and Fearing measuring ability. Such sensors are 
presently available in the form cf£ Radar altimeters{1*] and 
Laser ranging devices. Another possibility for a measuring 
method is a single sensor time sharing to cbhtain range and 
bearing tc kcth reflectors from a single device. Such a 


Sirgle senscr scheme is sketched in figure III-2. 


Once tke FWD and AFT parameters are available, they may 
then be used to determine other desired quantities. 
Subroutine HDGRAS (HeaDinG control for RAS) was develored to 
output tke desired heading corrected for heading difference 
@eeeship A and B and the prejected correcticn for distance 
errcr. This subroutine is listed in appendix A. The center 
Fange and tearing are the average of the FHD and AFT range 
and bearing cutput from Subroutine RBMEAS. 


The additional heading due to distance is projected as 


if ship E saintained its present course until it was 
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perpendicular to the center of ship A. The reasoning Eehind 
this is illustrated in figure III-3. If the present course 
will cause ship B to arrive on the station desired (DS), no 
heading chance 1S reguired. The expression for PSIACC (J 
Additioral heading due to Distance Correction) is: 


| pstapc = RSENS (DDC+DAeSIN (AA1) ) 


where: 
RSENS = Range SENSitivity gain 
DDC = Distance Desired Corrected for side of apprcach 
DA = center Distance Absolute (range) 


AA1 = 3€0 degrees - relative bearing of center position 


The beadince difference of ship A and B is desired 
Since, even if the range when alongside is correct, a large 
disparity in héading cannot be tolerated. It is realized 
that some heading difference (crabbing) is necessary to 
Maintain tke distance. This crabbing is due entirely to the 
pressure ferces modeled in chapter II. This heading 
difference is found by computing the difference in the 
perpendicular prcejection between the FWD and AFT 
measurements and finding the arcsin of this difference 
divided Ey the distance between sensors. Figure IIiI-4 


indicates a sample of this procedure. 


The expressicn for total desired heading 1s given as 
fEcllcws: 


PSIALCH+WNTSENSePSIDIFtPSIE 


PSIDES 
where: 


PSICDIF = 9 additional heading due to heading DIFference 


WISENS WeighTed heading difference SENSitivity qain 
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Heading Difference Calculation 
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PSIDES = § (heading) DESired 


Thretghcut the subroutines and main DSL programs, the 
Function CEGRAD (conversion of degrees to radians and 
radians tc degrees) is used freely. An explanation and 


listing cf this functicn are presented in aprendix A. 


“The angular velocity of the receiving ship's head is 
alec cf cencern in the RAS Situation. This quantity may be 
thought of as Similar to tachometer feedback in a Simple 
serve ccntrcel system; and is necessary to damp out the 
response (the responses of this control system without this 
feedback is presented in the latter section of this 


chapter). 


The desired rudder command is a combination cf the 
desired beading, angular velocity feedback, and a rudder 
gain as fcllcws: 


Desired Rudder = (YAWD2-PSIDEL+BDOTFSB) «GN 
where: 

YAWDZ = heading of ship #8 (in degrees) 

PSILED =ESIDES (in degrees) 

BDOTFB = VFBGeBLOTZD 

VFG = Velocity FeedBack Gain 


BLOTZD = angular velocity of ship B heading angle (in 
degrees/sec) 


BDOTFB = angular velocity FeedBack 


RGN = Rudder Gain 


The convention for rudder response dictates negative 
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rudder as kéing right rudder, which causes positive yaw. 


This necessitates making the desired rudder the negative of 
the forcing function and feedback quantities. The block 
diagram cf figure III-5 presents the control loop 
measurement inputs to desired rudder command. 


from 


72 





} 
! 
| 
{ 
| 
| 





| 
| 
81h 
45 oe | 
| | 
mene ke: 
= | rt "i | 
t | t 
és Saaisd\ < ae < | 
| + ft sf 
a ! | 
DadA 4Z ) p= 
| | 
\'/ ) _— 
| | 
Gurag | ae MS 
p | | 
-— ty 


ih 
530 | | 
P ia SYHOMH BINT LACYANS 


werbetq yooTq Torz3uU0D Tappny 
G-III eanbtga 





12 





Optimization 


Thus far the control chcice has identified four gains 
(RSENS, &ISENS, RGN, VFBG) that must be set for [reper 
position attainment. The necnlinear nature of the systen 
which afprears in the form of distance measurements, 
interactive forces and rudder modeling do not allicw for 
straight forward determination’ of these gains with normal 
Seersal cchtrcl theory. , 


Grossly nonlinear systems require special handling to 
determine prefer gain settings. The method chosen for this 
purpose was an optimization algorithm develored by M. J. Box 
(percgrammed locally as subroutine POXPLX). This subroutine, 
listed in appendix A, was used to locate the cost function 
saddle pcint in four dimensional space (the dimensicns being 
the previcusly mentioned gains). The drawback asscciated 
with this wethod is the necessity of iterating the complete 
ncnliinear sinulaticn within function FE for every evaluation 
of the ccst function. The gains sought were found, but 
unfortunately only after 2 1/2 to 3 nours of CPU time with 


every 400 iterations allowed. 


The weéechanics involved in optimizing the chosen cost 
function include required sub-calculations in many functions 
and sukroutines. Figure III-6 is a flow chart which 


demonstrates the steps, subroutines and functions required. 


Initial optimization was accomplished for one set of 
fiteradi conditions. By looking at the RAS situaticn, a 
rrckable set cf circumstances were envisioned. The scenario 


setting is the approach phase where the rerflenishing ship 
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starts a wide approach at 0.4 ship lengths (211 feet) 
lateral displacement and 5 ship lengths (2639.0 feet) astern 
of the suprly ship. The desired final position is alongside 
at 0.2.ship lengths (105.6 feet) lateral separation. Eoth 
Ships have tke same initial heading (YAW angle). The supply 
ship is at 15 kts. (1.0 normalized speed) and the receiving 
Ship makes its approach at 22.5 kts (speed control will be 
covered later in this chapter). 


Normai ccestinc cf displacement error with the integral 
of the squared error (ISE) was considered as the 
optimization tcol in subroutine BOXPLX. However, this type 
cf perfctmance measure would weigh the initial displacement 
errer equally with the final position error. This problem 
can be circumvented by comparing the displacement error to a 
pre-ccmputed reference track instead of to the desired 
displacement. For the envisioned scenario, it was conceived 
that the ccst function should weight the distance 
displacement teavier when the ships are alongside than when 


the apprcach is started 5 ship lengths astern. 


This wés accomplished by using the integral cf time 
times the aktsclute error (ITAE) as the optimization 
performance measure. The reference displacement was 
considered tke desired position displacement. The object 


furcticn can then Le written as: 


te 
CBJ -/ te|DD=-ADY] dt 


t, 
where: 
DD = Desired Distance 
ACY = Actual Displacement in the Y directicn 


t = time 
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A perfcrmance measure that is designed to obtain good 
performance must also take into account other factors 
besides just fosition accuracy. Consequently, another cost 
criterion was decided upon which would also set the gains to 
reduce tke agcunt cf rudder contrcl required when alongside. 
This particular feature is derived from the desire net to 
over contrcl with the rudder in such clcse proximity to 
ancther vessel. The inclusion of this term in the 
performance measure is weighted by unity while the distance 
accuracy iS weighted by a factor of 10.0. This will tend to 
allow rudder action if the desired pesition is not 
Maintained. The final approach phase cost function for 


oktaining cptimum gains has the form: 


t¢ 
CEJ = / te(10.Q0e«{DD-ADY|{[+1.Q0¢ej{D2]) dt 
te 
where the adéiticnal term is: 


D2 = rudcer response of the replenishing shir 


Cc. Resuits 


In the process of deciding on the best gain definitions 
previously mentioned, many optimization runs were made. 
Each set cf gains were then simulated in a corresponding DSL 
program to cktain performance confirmaticn. Many of these 
runs did not live up to expectations; causing re-evaluation 
of the ccntrcl scheme until the one presented in this thesis 


was formulated. 


Table III-1 shows the input upper and lower limits of 
search (EU, EL), starting value guess (XS), optimum gain 
settings (Output) and associated object function value (CBJ) 
for 20.0 seccnd normalized time simulation run in function 
ne These values were then introduced into the DSL 
Simulaticn pregram listed as program #5. The results of 
this sSsimulaticn are shown in figures III-7 thru IfI-12. The 
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Secaqraphic flct cf figure III-10 indicates excellent 
positioning in the lateral direction while the rudder 
response cf figure IITI-12 shows that it settles out to a 
fairly ccnstant steady state value as the ship settles into 
its desired fcsiticn. The time coordinates in all flots are 


shown in actual full scale time. 


d. Contrcel Testing 


2 





Now that the "proper" gain settings were cbtained, more 
extensive testing of the control system was required. Three 
different tests were contemplated: (1) allow a large 
perturbaticn turn of the reference ship (supply ship), (2) 
start approach of the receiving ship (ccntrol ship) from 
different initial conditions cf lateral and hcrizontal 
displacezents, and (3) induce external perturbations in the 


form of wave fcrces. 


The first test was Simulated by turning the reference 
ship by ncrmal rudder action of figure IiI-13. Tavs (ton 
With 5 degrees rudder accounted for a total reference yaw 
Change of 15 degrees. The rudder action of the ccntrolled 
ship shcwn in figure III-14 was as expected. However, the 
distance taintainment portrayed in figure III-15 was totally 
unacceptable. The maximum excursion frem the desired 
distance cf 105.56 feet (0.2 normalizes distance) was £5.419 
feet cc. 105 normalized distance). Variances of this 


Magnitude cannot be tolerated in the RAS environment. 


Faced with this situation, the tact chosen was to 
re-evaluate tke gains for the new scenario which is called 
the turn phase. In this phase the initial conditions assume 
steady state pesitioning alongside such that the lateral 
fpesition displacement (DY) is equal to the desired distance 
(105.56 feet (0.2 normalized) ] and that the hcrizontal 


Ecsition displacement (DX) is O20 (alongside). 
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Scre initial ferturbation is introduced by assuminc the 
relative yaw angle when alongside is negligible. 


The same procedure was followed in obtaining gains that 
would optimize a chosen cost function. Figure III-6 still 
applies excerpt that function FEA is replaced by functicn FEB 
{listed in aprendix A) to simulate the new conditions. 


Cost function criteria change in this instance since 
the ships start at the desired pcesition and optimally stay 
at the same relative positions. Also, the rudder response 
to such a large turning perturbation must be free to cause 
achievement Ci the desired position. Due to these 
considerations, the integral of the absolute error (IAB) 
performance measure was chosen for the optimization 


criterion ané can Ee written as: 


ty 
OEJ = / {ADY|{ dt 
t 


o 


Table III-~z shows the results of the turn phase 
optimizaticn and the comparison with the approach phase 
gains. Acain CSL simulation was performed using the turn 
phase scénario. Pigures III-16 thru III-21 portray the 
graphical results. The rudder response cf figure IIfI-21 
indicates very sensitive response to the interactive forces 
shewn in figtres III-17 and IfI-18. The lateral distance 
separaticn cf figure III-20 indicates excellent fosition 
Maintainment with maximum excursion error of only 2 feet 
(0.0038 norzalized). This minimal variaticn is well within 


that which can be tolerated in the RAS envircnment. 
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f. Centirved Control Testing 
a Per a ap 


—~— Se a Se a a == 


To alleviate suspicions that the response frog the 
gains oktained in the approach phase could be improved by 
the gains ckEtained in the turn phase, a simulation of the 
approach frhase was accomplished with the new gains. Figure 
III-22 is tke graphic display of the effect of these gains 


on the approach phase lateral distance positioning. 


Careful analysis of the results thus far clearly 
indicate the need for an adaptive control scheme-to allow 
gain adafptaticn to meet the design specifications. A _ full 
acaptive ccntrol scheme for systems of this type is outside 
the scope cf this thesis. References 15 thru 22 are 
‘indicatiors cf some of the literature available for ptrsuit 


of a ccmpletely adaptive control systen. 


What was done here is develcpment of a Simple algcrithao 
to sense when the cenditions were adequate to switch fron 
cneé set cf gains tc another. This may be done with the two 
sets of gains developed thus far. However, for the sake of 
Sigulaticn efficiency, a third set of gains was introduced. 
This third set amounts to a change of one approach gain 
(RSENS) which has previously been defined as the range 
sensitivity cain. The simulation efficiency is increased by 
decreasing tke time required fcr the approach phase to reach 
Steady state. A conseguence of this procedure is a 
Teenforcement of the need for a completely adaptive ccentrol 
scheme. 


Repeated simulation revealed that ccmmencing the turn 
(in effect switching gains), before a reasonable steady 
State was reached caused results similar te those shcwn in 
figure IIIJ-zzZ. An increase in RSENS to a value of 4.0 when 


the lateral separation error is less than 0.05 (normalized) 
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and greater than 0.005 (normalized) forces acceptable steady 
state in afrroxizmately 1/2 the time previously reguired 


usirg a single set of approach gains. 


Subroutine S&TCH (listed in appendix A) inccrpcrates 
this simple acaptive gain schedule with a ccunter tmechanisn 
to sense when steady state 1s reached. Further study 
indicated a need to damp the yaw oscillaticns to a greater 
extent if the yaw velocity (BLOT2D) exceeded 2.0 degrees/sec 
when the cairs are initially switched to the more sensitive 
enes of the turn phase. This is an artifical adaptive gain 
fer VFBG caused by computer time restricticns prevalent in a 
full scale computer simulation where both the approach and 
turn phases are desired. If the gain switching point is 
moved up in time, as would normally be the case in a real 


life situaticn, this damping increase would not be required. 


The results of the full scale Simulation using cogputer 
program #6 are shown in figures III-23 thru IIIrI-34. The 
aprroach fhase flots. of figures III-23 thru ITI-2@ show 
definite improvement over that previously shown in figures 
at / ereu ITii—i2. Figure III-27 indicates that the 
Oovershcot is reduced to 10.6 feet (0.02 normalized) as 
opposed to 113.9 feet (0.034 ncecrmalized) that was prevalent 
in figure III-11. 


The ttrn phase plots are shown in figures III-29 thru 
fitI-34 and shcw responses very similar to those shcewn 
previously in figures III-16 thru III-21. The only 
Significart differences occur in the initial responses which 
are due to the incorrect initialization when the turn phase 
was Simulated individually. 
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Men Varying initial Conditions 


The results oktained in the previous section are most 
gratifying fut actually incomplete. This system must work 
fer cther initial. conditions quite different frem those 
envisioned in the optimization scenario. The initial 
approach can realistically commence at points other tkan 5 
ship lengths astern and displaced by 0.4 ship lengths. 


By sSigrulating this system with varying initial 
Ecositions, the relative efficiency and wortk of the ccntrol 
system can te cbserved. This was done in successive test 
runs whose iritial conditions and corresponding plot figures 
are tabkulatec in table III-3. For the sake cf brevity only 
these figures required to illustrate the relative efficiency 
of the ccntrcl system are included. The correspcnding 
initial cptinizaticn simulation figures are listed for cross 
reference. The turn phase plots for all runs except 4 and 6 
exactly match that of the initial simulation and are not 


repeated kere. 


Runs 3, 5 énd 6 were accomplished to show that no 
ambiguities exist in the control scheme to prohibit adequate 
real life initial conditions. Run 3 simulates the situation 
mcst often encountered by this author in the RAS 
environment. This scenario starts the ccecntrol shif dead 
astern at 5.0 ship lengths and brings it alecngside at 0.2 
Shig lengths lateral separation. 


Run 5 1S a situation where the apprcaching ship is 
purposely placed out of position on the wrong side for 
approach. The control scheme adequately corrects the 
placement error and will do so for all cases of this type, 
Prcevided that there is adequate maneuvering room astern of 


the reference vessel (in this case 2.6 ship lengths was 
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experienced which gives 1.6 ship lengths bow to stern 


clearance). 


The purrpese cf run 6 is to frovide simulation fcr an 
approach from the opposite side again disproving any ccncern 
for ambiguity in the trignometric measurement scheme 
utilized. In all runs it must be emphasized that DD is the 
positive absclute distance desired and that IS provides the 
code flag for the desired side of approach. The system will 
work with DD set to some negative guantity; but the side of 
approach *ili reverse itself and the positicn placement will 


ke correct, Eut cn the side not desired. 


Run 4 takes the desired distance in to 0.15 ship 
lengths (&€0.0 feet). This distance is usually the minimun 
desired by a prudent seaman. Again, even with this minimun 
distance, the control system performs uf to desired 
Standards. Tke importance of this. run cannct be overlcoked. 
Performance cf the system at this extremum indicates that 
the gains utilized are correct for all expected conditions 
encountered in calm _ seas. Figures III-59 thru III-64 
pertray tke remaining plots obtained in run 4. 


Icrmance in Sea 


h. Perfcrmance in tate 


The calm sea performance of the heading control systen 
is only part of the system testing required. Of even 
greater ccncern is the adequacy of the control when sea 
state 1s intreduced. Section D.2. of chapter II models the 
three ccuponents of waves with two sinuscids and a small 
random impulse wave. These forces were intrcduced intc the 
tctal RAS simulation as shown in computer program #7. [In 
this program the wave length (WL) is set to cne ship iength 
and the wave direction (WD) is -015 degrees true. This 
scenario allcws for a port turn into the prevailing sé€a as 


1s commen practice in experienced RAS evolutions. BY 
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Binimizing the perturbation forces on yaw and lateral 
direction, a smocther RAS can be accomplished thus aiding 
safety and ccmfort during the actual transfer. The wave 
force maximum is taken as 0.05685. Runs were simulated 
“which used maximum wave forces in the range 0.1137 to 
0.65685, wave lengths from 0.5 to 1.5 ship lengths and wave 
directions 015 to -015 degrees off the initial replenishment 
course. The ccntrol system handled all of the perturkations 
well except fcr the cases of a wave length of 1.5. This 
length of wave with a. force of 0.05685 exceeded the ccntrol 
systems capability in that the steady state conditions were 
not met before a turn was commenced. Figure III-65 shows 
this instability in the lateral distance DY of the turn 
phase. It is felt that the modeling inadequacies of the sea 
state develcgyment of chapter II coupled with a simple 
adaptive gain scheme are the source of the problem. This 
Same phercrenon is covered in greater detail in the 
longitudinal fosition offset testing portion of the velocity 
control section of this chapter. 


Proctlems of this type also manifest themselves in some 
cases when the wave force maximum (WFMA) was close tc the 
0.1137 value. If the sea state beccmes excessive, which 
this value represents, a different gain schedule or, at 


best, a mere complex adaptive gain scheme is called for. 


The flets produced by computer program #7 are presented 
aS a representative indication of the effectiveness cf the 
ccntrol system in the presence of a sea _ state. Figure 
III+66 gives the yaw results of the appreach phase which 
indicates the effect of the wave action. The correspcnding 
rudder acticn cf figure ILI-67 compensates to give the 
smooth lateral distance shown in figure III-68. The wave 
prcfile is Shown in figure III-69 with curve 2 being WY and 
curve 3 -éinc WN. Curve 1 is the WX profile which was not 


used in this run but will be utilized in the speed ccntrol 
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section later in this chapter. Similar curves are portrayed 
for the turn fhase. Figure III-70 is the yaw difference 
between the two ships (remembering that the reference ship 
is nct being perturbed by the interaction forces or the wave 
forces). Figure III-71 is the lateral distance DY 
maintained Ey the rudder response of figure III-72. The 
maximum lateral separation in the turn phase is 0.0037 ship 
lengths (1.S5 feet). The wave profile is shown in figure 
III-73 with the same wave force curve sequence as the 


approach rfhase. 


As can ke seen from these flots, the ccntrol systen 
Operates very effectively in the presence of a sea state. 
Again, the develcfment of a much more complex adaptive gain 
scheme is required to allow exceptionally high sea state. 
It is felt that the control system presented in this thesis 
is adequate for most Situaticns that are encountered in the 
RAS envircnment. Only the extreme perturbations that chance 
would allow must be accounted for in a more complex adaptive 


gain scheme. 
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B. VELOCITY CONTROL 


One advantage derived from using the linearized 
equations cf motion is the decoupling of the velccity 
compenents frem the remaining equations cf motion. This. 
allows separation of the design procedures for Lateral 
separaticn ccntrol and velocity ccntrol. Section A of this 
charter designed the Lateral separation ccntrol using the 
Simple speed contrcl algorithm shown in figure IITI-74. This 
control cutput was used directly as the ship's speed (CDCT2) 
in the mcdel simulation where no attempt was made to use the 
engine respense developed in chapter IIT. Function SPDCTE of 
appendix A shcws the control used. 
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Figure III-74 
Non-optimum Speed Law 


Because of this decoupling assumpticn, any valid 


approach speed ccntrol can be used, if used consistantly, 
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for such a design. However, in the RAS environment, 
cctplete disasscciation is not possible. Recombination 
occurs in the interactive forces and moments which depend 
upcn the longitudinal distance as well as the lateral 
distance, Consequently, speed, Which is directly 
responsikle fer the longitudinal distance, has a direct 
relation to the lateral distance attainment and 
maintaingent. 

The remaining farts of this chapter deal with the 
development cf a viakle speed ccntrol algcrithm and the 
testing cf the designed systen. 





i. Type cf Centrel 


Whenever two ships maneuver for replenishment at sea 
(RAS), the prime ccnsiderations are the time required for 
approach and the accuracy of position keeping flus 


conservaticn of fuel. 


The ncnlinear control law of figure III-75 is designed 
to maintain a preselected approach speed for nininup 
approach time. The proper location of the switching point 
increases the complexity of the soluticn since the tise of 
Switching frcem this speed is determined by the dynamics of 
the nenlinear fosition attainment loop, Once this pcsition 
1s reached, the speed controller is switched dawn toa 
linear perticn of the control law to allow control for 
perturbations about the operating position. However, small 
perturbaticns about this operating point can be tolerated 
and, in fact, are desired to allew for conservation of fuel. 
Selection of this dead zone is wholly dependent on the 
accuracy required LJCeEInal position. Figure f{II-75 
indicates a dead zone extending to +0.001 nermalized 
distance which in this case translates to +0.53 feet. 


Systems fcr which fuel considerations are net a motivating 
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factor say te designed without this part of the contrcl law 
ta allow finer tracking in the position loop. 
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The speed control law as explained abcve is shcwn in 
figure III-75 for an initial approach speed cf SPDOQ2 and a 
final estimated reference speed of SPDO1, with ADX being the 
dynamic fosition feedback defined as the longitudinal 
distance between centers of the ships referenced to the 
controlled ship's heading. Analytically, the linear pcrtion 
of the contrcel law is written as> 


SEDCTR = -ADXe(SPDC2-SEDO1) +SPDOT 


Symetric continuation of the ccntrol law accounts for 
Operation cn both sides of the operating point. 
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2. Cptimization 


Using this euch simplified model of chapter IIT and the 
basic ccatrel law of figure III-75, the desired switching 
curve can be established. An optimization subroutine such 
as Subreutine BOXPLX can be used to iteratively oktain the 


optimum switching rosition (Sd). for representative initial 


approach speeds.u, Figure III-76 is a flow chart cf the 
subroutines and functions required for speed ecntEol 
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Figure III-76 
Optimization Flow Chart 


Optimization. The major merit of this nenlinear contrel law 
Stems frcz the predetermination of the switching pcint for 


all possible conditions of initial speeds. This arfriori 











or ee rir ee 








knowlege allows for offline computation cf the switching 
pesiticn prior to commencing the approach. The cost 
function used for Se timizac.on eS the DEAE which 
accomplishes two objectives. First, it forces the appfroach 
tc be acccmplished in minimum time. Secondly, it insures 
that the fuel expenditure will be optimized in the 
elizminaticn cf most overshoot and bang-bang control in the 
dead Zone fortion of the control law. The final value of 
the poSiticn error must be within the specified dead zone 
and the ternzinal speed must match the reference speed 
(SPDO1). The cost function has the following form: 


ty 
s- f (te JADX|) dt 


+, 
Table III-4 is a comparison of the optimization fcuns 
With various initial speeds. The values shown for S& must 
be multiplied by the speed differential (SPDO2-SPLBO1t) to 
obtain the ccrresponding value of ADX. The max/min values 
show the kand cf£ values which preduce the optimum cost. 
This cance cf values is attributed to the integration ster 
size used in the optimization program. Experience with this 
particular optimization program indicates that errcneous 
values of the switching point are found if the step size is 
nct carefully chosen. The step size may be adequate for 
Macegraticn, but net for location of the switching foint. 


The fpeints obtained from the optimization runs are 
Fletted in figure IiI-77. These points define the nenlinear 
Switching curve which must be stored in the computer to 
insure optimal operation of the speed ccentrol for all 
approach speeds. From here there are many frocedure oftions 
Ofen. These cptions have as a goal some usable fer for 
predicting the optimal switching point for any set of 
initial ccnditions. One may choose linear straight line 
segments with an interpolation routine, cr a closed forn 


Switching curve pclynomial. Due to the availability of a 


160 








INITIAL CURVE POINTS 
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Table III-4 


Optimization Results 
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hykrid ccnfigured XDS 9300 digital computer and AGT-10 
graphics terninal, the latter course was taken. 


A pclyncmial curve fitting algorithm was used to cbtain 
the reguired pclynomial coefficients of best fit. This was 
done for fpelynomials of order 1 thru 5. The coefficients 
and the sum cf the sguares of deviation frem the original 
Feints are tabulated in table III-5. The selection cf the 
order to te used is highly defendent on the degree of 
accuracy required. In the RAS problem, the average error 
introduced for a first order fit is 8.0 feet(1.07? séc), 
while the fifth order fit introduces an average error of 
1.35 feet(C.180 sec). Prior acceptance of errors intrcduced 
by an integraticn (and problem) step size cf 0.8 sec allows 
for use cf a second order fit without any degradaticn of 
Simulaticn accuracy [second order average error is 2.848 
feet (0.38 sec) J]. The graphic display of figure IIiiI-78 
indicates very little difference in the switching curves for 
second tc fifth crder polynomial fits. For the sake of 
accuracy, and owing to the computer control methods of this 
thesis, tke fifth order polynomial fit shown separately in 
figure III-79 is used for determination of the switching 


Foint locaticn. 





A trte test of the control law 1s accomplished when it 
is introduced in a computer program fcr a complete RAS 
Simulaticn. Considering the performance of this ccntroller 
in a complex envircnment of full scale RAS simulaticn allows 


maximum verification of the controller design. 


The scenario for this simulation initially fositions 
the ships such that the ship being controlled starts an 
approach 5 ship lengths (2639 feet) behind the reference 
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ship and displaced 0.4 ship lengths (211 feet) to the right. 
The desired final position is alongside and displaced 0.2 
ship lengths (106 feet). The heading contrel system used is 
developed in section A of this chapter. 


The approach phase 1s accomplished with the _ speed 
desired and speed acquired shown in figure III-80 with the 
corresponding position attainment exhibited sh ial figure 
ier 8 1. These plots show excellent switching and oftimal 


position attainment. 


The next step is to insure that the fosition keeping 
lecp will maintain the desired position with an induced 
perturbation. This is accomplished by turning the reference 
ship away frem the control ship a total of 15 degrees to 
observe the reaction of the speed control locp. The 
reference ship's turn causes the relative mction between the 
ships to ce altered, making the control ship lag the desired 
position. The nonlinear control system is designed to 
correct this situation as soon as the actual position is 
outside the limits of the dead zone. Figure III-82 displays 
the desired speed and acquired speed for the control ship. 
Pigure III-83 indicates that the corresfpcnding fosition 
deviates frem the desired by 0.0154 ship lengths (8.13 feet) 
at the maximum excursion. This is weil within the limits of 


acceptakrility for such a drastic perturbaticn. 


The introduction of velocity control was acccmplished 
by combining the simplified engine response of chapter II 
and the speed ccntrol law developed here. By setting the 
speed desired (SPDDES) equai to the output of Function 
SPEREC and scaling the speed error (SPLERR) to the 
nondimensional equations of motion, the velocity lcop is 
initiated. The auxillary equations added to those presented 
in chapter II are: 
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SPTDES = SPDREC (ADX,SPDO1,SPD02, Sh) 
CDCT2 = INTGRL(U02,SPDERR*¥LUC) 


These equaticns are introduced in computer program #8 to 
produce figures I1I-80 thru I1Ii-83. 


Further system study indicates that the reference ship 
speed must be known to a fairly high degree of accuracy. 
Without afpriori knewlege of the reference ship speed, a 
constant bias is introduced. The amount of bias allcwable 
defines tke fermissible uncertainty in the reference ship's 
initial epfeed. This bias can amount to as much as 0.1 Ship 
lengths (€4.48 feet) for a reference speed inaccuracy cf 2.5 
kncets (0.1 normalized speed). However, it is felt that the 
reference speed in any practical situation will initially be 
kncewn to within 0.5 knots (0.02 normalized speed). This 


mcre practical error will introduce a bias of only 16 feet. 


Cther feedback parameters can be used to offset the 
lack of apricri kncwlege of the reference ship speed. Since 
the reference ship is tracked with a high accuracy rance and 
bearing dévice and the controlled ship's speed is seasured, 
a decoupled multivariable scheme is used tc further refine 
the reference ship speed. With high resciution devices 
presently available{[t*], 1t 1s estimated that this can be 
dcne practically to within 0.05 knots (0.002 normalized 
speed). This would bring the offset bias to 1.6 feet; well 
Within previously defined errors introduced by integration 
step size. 
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Thrceughcut the development of the heading contrcl and 
speed contrcl, the scenario has followed the conditicn that 
tke final position would be longitudinally alongside. 
Although this is a good assumption for ships of the _ same 
type, it dees not account for RAS station differences for 
different ship types. To alleviate this disparity, function 
SPDREC was redesigned to allow pre-planned cffset condition 
to exist. Function SPDOFC cf appendix A is a result of this 


redesign. 


Simulaticn runs, with a change of the speed ccntrol 
function only, resulted in some unstable ccnditions existing 
in the headirg ccentrol loop. The cause of this pheéencmenon 
stems back to the adaptive gain scheme used and the changes 
made te fcrce the control loop te a steady state value prior 
to aturn. Ey using a favorite ploy of experienced ccnning 
officers, this problem is alleviated. The ploy is to take 
the ship alcngside and then either drop back to station or 
surge forward to station. This method is accomplished by 
setting the initial offset (XOFS) to 0.0. The final desired 
offset (XCFSD) is stored and not used until the ship is 
settled out alongside. It is subsequently used as shcwn in 
the follcwing Fortran code: 


IF (ATIME.GT.450.0) XOFS = XOFSD 


This methcd solved the gain transition problem. It did 
not, however, give a completely stable simulation Cun. 
Unstaktle ccnditions still existed at the end of the turn 
phase. This is not suprising, considering the heading 
ccntrcl optimization method used. The set of gains 


previously fcund were for the alongside scenario only. 





Different interactive forces and moments at the cffset 
position cause these gains to be no longer cptinal. 


By relaxing the control loop in the heading velocity 
feedback gain (VFBG), sub-optinal control at all practical 
orfset pcsiticns 1s achieved.j The gain VFBG was changed 
frem (0.064028 to O.1 in the turm phase adaptive gain 
schedule without significant less cf centrol efficiency for 
alcengside operaticn (2.3 feet maximum excursion vice 2.0 
feet previously obtained). Subroutine SWICHF of appendix A. 
reflects the gain change and offset calculations required. 
.Ccmputer program #9 incorporates the changes required for 
offset simulation. Table III-6 is a cross reference Listing 
of the plets obtained. From these figures, the effect of 
different Icngitudinal positicns is readily apparent. An 
offset cf 0.1, eguating to 52.8 feet, causes greater Lateral 
excursions when astern (XOFSD = -0.1) of the alongside 
position than when ahead (KOFSD = 0.1). The Longitudinal 
Bosition naintainment, however, is essentially the same in 
all cases. | 


Se 


nor SD. 


Lateral 

Distance DY 

Yaw 

Difference 

Speed 

Response 

Longitudinal 

See ele lete tate 


: 
Note : These plots are the same as those obtained from 
computer program #8 and are not repeated here. 
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Approach Phase Plots Turn Phase Plots 
A - B e A B e 








Table III-6 
Position Offset Testing Cross Reference 
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An alternative to the —peitat shown here 1S again a 
completely acaptive gain scheme which would achieve oftimal 
control instead of the sub-optimal control settled for here. 
The alternative may become even more important if the 
ncnlinear terms of the equations of moticn are considered. 
This would ccuple the heading and speed control designs to a 
larger extent than encountered in the interactive forces and 


gcrents. 
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The final testing procedure involves validation of - the 
speed ccntrcl system in the presence of waves. This 
perturbaticn testing continues that started in section A. of 
this chapter for heading ccntrol. In chapter II the WX 
force was acdeled thru the intermediate force IF32 as: 


* 


IFiZ = KC1eD2+NC2+KC1¢HX 


By intreducing the force in this way, a severe limitaticn is 
placed cn tke magnitude of the force. In the mariner model 
used, the KC1 coefficient (XDELR) is considered negligable 
cr, at best, cnly 0.00005. This translates, in the original 
equations cf motion, to a maximum speed perturbation of only 
O0.C255 kts. for the wave amplitude chosen. The second 
drawkack cf this method, with even greater ccnsequences, 1s 
that the perturbation is introduced before the contrcl loop. 
Delay of the wave perturbation is produced making it out of 
phase with the cther wave force {WY) and moment (WN). 


eee ae ee a ie ee _— — —_ oo — — _— _— — — 
- « -_— - - —_ s -—— -- _ -_-- _ — - - 


SPDOES 





Rigure: Eit=-102 
Block Diagram of Wave Introduction in Speed Loop 


MmimGLdeE to Ding about uniform introduction cf this 


wave force, its effect is inserted just past the integrator 
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of the speed contrcl loop as shown in figure III-102. This 
is coded in the DSL simulation program as: 


CDOT2 = INTGRL(U02,SPDERR*LUC) +KS2*WX 


A value of ~-1.0 for KS2 will give a maximum wave 
perturbation of 0.85275 kts. (a auch tore realistic 
perturbation for the high sea state simulated). Figure 
ItTI-103 portrays the speed desired and speed acquired for 
the apprcach phase in the presence of sea state. From this 
it can be seen that the speed acquired is very dependent 
upon the sea state present. The contrel law, however, 
presents a very stable reference for the speed locp which 
gives an approach longitudinal position (DX) plot 
indistinguishakle from that of figure III-81. More 
preminent perturbaticn results are evident in the turn phase 
flets cf figctres III-104 and IIi-105. The speed response of 
figure III-104 allows a maximum longitudinal position 
excursion of 9.5 feet (0.018 normalized position). as 


ccBpared with 8.28€ feet (0.0157 normalized) in calm séa. 
These results show that the speed ccentrol system is 


very stable and corrects well for large external 
perturbaticns. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLEUSICNS 


The results of this design study have been most 
gratifyirg. The kasic concepts initially ferceived for the 
RAS contrcl have been realized. The decoupled ship ccntrol 
in the RAS environment is a viable and plausible idea. This 
thesis ccntains a workable system for implementaticn of 
ccmuputer ccntrolled RAS. The achievement of 2.3 foot 
maximum excursion for lateral distance while both ships are 
in a tturn, and longitudinally offset by 53 feet isa 
phenomenal achievement. Having this kind of accuracy in RAS 
Ccperations, can vastly increase the safety of this 


ccaplicated and dangerous maneuver. 


The appreach phase of RAS can be a very hair raising 
experience. Night replenishment and sea state ccmplicate 
the "seargan's eve" method now employed in the fleet. Having 
a system that autcmatically handles the apprcach regardless 
of the adverse ccnditions can, again, de nothing more than 


increase the safety of the RAS maneuver. 


Schemes for ccmputer control of nonlinear systems and 
the purpcseful introduction of nonlinear ccntrol laws are 
beccming mere practical with the technological advances in 
Micro precessors. The ever increaSing number of U.S. Navy 
ships ielee in computer systems installed, makes digital 
computer shir control realizaktle in the present time frame. 
A gcod wicrce ccmputer or an existing installed computer 
(Such as cne used for the NTDS system) can be used in this 
veln. Frocurement of the hardware required for this RAS 
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system can be dissipated over time periods contingent cn the 
funding available. The supply ship requires cnly two 
reflectors fcr the range and bearing devices stationed on 
the recéiving ship. All ships can be outfitted with such 
reflectors at a minimal cost, while the bulk of the hardware 
can be introduced to the ships at regularly scheduled yard 
periods. 


In tke initial conception of this thesis, a secticn was 
rlanned fcr cpen ocean maneuvering. After scme research on 
this facet: of ship control, it was determined that werk in 
this area has already been documented[ ** ]f 25 ]. The 
existence of NTDS outputs for station attainment and single 


ship contrcl systems, made design in this area redundant. 


The concept cf integrated centralized ship contrcel has 
been in the tackgrcund for over a decadef 26]. Although 
given a lcw pricrity due to funding considerations, its 
implementaticn seems to be just arcund the corner 27}. 
However, a review of ref. 27 indicates a lack cf RAS 
capability. Whether this is an oversight in the article or 
neglected in the désign criteria is unknown. If it has been 
neglected in the design, a very real prceblem has fceen 
overlooked. The recent incidences of ship collisions while 
conducting RAS [28] emphasizes the need for inclusicn of 
this very cangerousS maneuver in the “Integrated Eridge 
System." Lack of technology can no longer be used as an 
excuse, This thesis and other research reports [291] have 
advanced the irzplementation feasibility tc a level that 
cannot be ignored. With these projects finalized in 
practical terms, their incorporation into fleet use is the 


next imperative step. 
A major effort in this area must be made. The ever 


increasing ccmplexities of today's naval shifs and the leads 


being placed cn the officers and men are such that ccmfuter 
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ccentrol must be used; and used now! We cannot afford the 
luxury of time to prove these system u crth, but must make 
ccncerted effcrts to get them implemented before the lives 
of 300+ men are lost. 


Whenever a complicated system such as a ship in the RAS 
situaticrn is enccuntered, many facets have ve be 
concurrently analyzed. This fact has caused inclusion of 
many diagrams in this thesis to illustrate the total 
picture. Each run, with a different ccndition, requires 
Many picts tc analyze the differences in the responses and 
the causes cf the differences. The computer programs shown 
de not reflect the actual run times in the JCL shcewn. As 
Many aS twenty flots were output in these programs in the 
tines listed. Analysis of the actual computation times show 
that the algcrithms run considerably under the time required 
for real time operation. The sampled data rate used in the 
Simulaticns was 0.11 seconds. This is well within the 
realizable data rates available in even the slcwest of 
today's ccmputers and microprocessors. The thrust of this 
consideraticn is that there are nc prceblems envisioned in 


ccnverting RAS simulation to real world RAS control. 


BP. RECCMMENETATICNS 


In the heading control design section of this thesis, 
the need for a completely adaptive gain scheme was cited. 
Again in the velccity control section, when a longitudinal 
offset was introduced, this need became even more evident. 
The first and most important recommendation for further 
study is the development of just such an adaptive gain 


scheme. 


The linear equations of motion should te replaced with 
nonlinear equations to validate the control designs advanced 
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in this thesis. Along with this, the hydredynamic 
coefficients for the Navy's modern ships are required to be 


able to design these control systems for todays vessels. 


It 1s further recommended that a ccncerted effcrt be 
made to cktain data on the interactive forces and mcments 
between shirs of dissimilar tyres and sizes. These forces 
and moments gust also be available for sea state conditions. 
In fact, tke whcle area on séa state effects cn the various 
ship types in the RAS situation and in open ccean 
Maneuvering needs attention. Not enough data was available 
fer this researcher to be able to pinpoint sea state effects 
on ship hulls. Since replenishment at sea is rarely 
ccnducted in the sterile condition of calm sea, these 
consideraticns are of utmost importance to allow testing of 


any contrcl system in the Simulation stage of development. 


The intent here is not to imply that the ccntrol 
systems pcertrayed in this thesis are the best for the RAS 
scenaric, but that the procedures used can te applied to any 
control scheme desired and bench marked to the ones 
ccrtained here. AS previously mentioned, much meaningful 
research and éesign must be accomplished tc allow system 
reliability and, mcre important, system acceptability Ly the 
officers and men who will ultimately trust their lives to 
it. This is a task that must not be taken lightly. 
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APPENDIX A 


Due to the lengthy nature of the computer prcgrams 
presented in this thesis, many functions and subroutines 
were develcpred tc simplify their presentation. This 
appendix lists these functions and * subroutines in 
alphabetical crder. The computer programs reference this 
appendix and indicate the placement of the required 


functions and subreutines. 


A Erief description for each listing is given to aid 
the reader in determining their purpose and use. The 
following is a listing of the functions and subroutines 


contained in this appendix in the order presented: 


SUERCCTINE BOXPLX 
FUNCTION DEGRAD 

FUNCTION DELAY 

FUNCTICN FE - RUN A (FEA) 
FUNCTICN FE - RUN E (FEB) 
FUNCTICN FE - RUN C (FEC) 
SCERCUTINE HDGRAS 
FUNCTION KE 

MAIN PROGRAM FOR FUNCTION MINIMIZATIONS (MINIEXPX) 
SUEROUTINE RBMEAS 
FUNCTION RKLDEQ 
SUEROUTINE SLOPES 
FUNCTION SPINIT 

FUNCTION SPDCTR 

PUNCIICN SBDOFC 

FUNCTION SPDREC 

FUNCTION SWCL 

SUERCUTINE SWITCH 





SUERCUTINE SWTCHF 
SGEBCUTINE TRANS 
FUNCTION XLIMIT 
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SUBROUTINE EOXPLX 


This sukrcutine was used fer all optimization runs in 
heading ccntrcel and speed control. It was programmed 
lecally ard is part of the IBM 360 SSP library at the Naval 
Postgraduate Schocl. A full explaination and description is 
shown in the first few pages of the subroutine listing. 
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SUBROUTINE BOXPLX 


elelelalelelelelel ele) elo) e*loleloelolelelelel®lolole)2) el elelelele\ele)e]elol@ el elele) 
MI ALT SUD OE OD OMA AT LOM ODOM OUM FIND MONO BAIMFIAOM ODO SHAM 
ODDO OOOO OC ee eet es eed cd peed eet ed eet et CIC QU CUI CUCU AIR AICS) CO CAH MAN CAO OS SO 
ODODDDOOWO OOOO QM OO OO DODO OOCGCODAOD OOOO COOOCOCOOCO9CO000 
>< >€ OK OE OK Ok OK OK OK OK OK OK OK OK OK OK OK OK OK OK’ OK OK’ OK DOK’ OK’ OK’ OK’ OK OK’ OK Oo’ OK’ Oo >’ Oo’ OK’ OOK’ OC’ OOo’ OO 
Qa2aQaaandcaaagaadcaagadaadaaaadgaaaaaaaadaaaadaaaaoaa 
>< >< >< Oo OK Oe Oe OK OK OK’ OK OK OK OK OK OK OK OK’ DK’ OK OK Oo” OK OK ot OK Oo’ OK Oo’ OC’ >t Oo’ O>*K’ D>’ OK’ OK’ O>*K > OOK’ Oo’ OK OS. 
MOMMdDODODOMDMMMD ODD ODDDADOHDODDODDDDADMODDHODDHODAHODHOADODHOHO® 


®eeeoeeaogeesoaoeaeaoaneeevep@e@eseeeeseweaeeeoeseowsceeeseeeoaseseeoaeegweoeaeneo~nseeosae@aeo@e@swseoeao ee eoeeedceane © 6 


HO) 


(CATEGORY 


SUBROUTINE BOXPLX 


PURPOSE 


Ta) 

Ld 
Oo > 
Z2= @M mi Ww Te) 
ead Iic> Mw 
bo FNL ee << 
LOK Xyv>RKHWTa 


Quel CWO <n 
Orem SeWW O>a 
WAWIWM GMM wW 
Sire RH Om Ww 
wLme2Sz OD STW 
OWA tl 2 bk 
Z2atooqow «— 2 
SWwOC]Z TZ eww 
UW > Ore WU) bk > 0 oY 
w= WAZ WI<Wwd 
CO fae fee 7 LO WL OT ke 
OOOO > Os 
CWOOQOZe tm Ee 
a3 ae | wll 
Muawkm Zid De 
WO a=W eft oe ZO 
XK TPUTLPSA zee 
eet NOE 
fe a me WANS 
wtrvQ2arvOoOrwvw<MN 
=o Ox§€OTe Mew 
wheel =o - Or 
Om OG N—eNMNWO 
WOZ Oe Ste Seu 
wt e2oamtZaZzaac 
Ont «Deez ff 
ee POOAYQSSI ae 
ZJmaQOnmo Orrxwnmn 
O<arW LL O00 be LL COLL 
Um Or MOM Wi 
ALS tw OZ Wwi-—cw 
POa e2>DOWZ—t 
MEOKTN OW) te 
Ud mili Jae we «=O 
as kh ae 2< 
=— ~§$OQONO=erae> 
eae Zee Ow 
Sait O ef SS pee YALL 
OxxvO sOASIDYZ 
Cte TeyvTtva wl 
Ozawa DS wad 
> OQwakeWwe=T 2 
Not ~~ Eawn= hm WwW 
oat DBZO00Q —aow 
Iw « OBSSZ> W 
OuwnOre<f azo 
DS ew ad <f Oz 
tong Qa. ny ZO 
ZEKE STO— eOW- 
>< ee) C3 LL SD Yt ey CLL 
wa ZW tu JOL 
Qa we SO 
MS MOWLA@mZaw 
Sy > ewe at S| DLL 
O<f Nike Ct es LL I 


CALL BOXPLX (NV NAV NPRa NTA ReXSpTP2XUsXbeYMN, IER) 


USAGE 


DESCRIPTION OF PARAMETERS 


ce <crM 
LL) <E fee fa fae 
(ius = 
oe es eee 
—Oh me SP 
zai tu 
vod re LG fe 
Tree ei )<f 
ZtOWW 
<I> ZNawM 


> 
ze 


> Zl - @ 
Ud ZNA-™ 
~> TL pws 
YSZoOmO NO 
{Oct M2 J 
mit) <i LL 
cael °JZUud 
ved LL LL CLL tt beg 
M=Z=OSIK< OC 
={O a0w et 
qt We den2Z 
N= — AU me 
ee a eee 
oretoomw 
>> Wao 
fae a Lh) beet eee 
WO LL > Wc 
Owrrsanatw 
= -tRKO>Y 
SU m—<tr 
ae (Ud LL ee Se 
merit oO 
Lt oe tL <t 
Tr hd UL, > bee CS) Pt pee 
meOWw<Et = JO 
CQ AaKnKnK~aZ 
OO 2Zewx<D> 
zeOataD DD 
— e Oiuw<f 
ZN gf 2 
m= o<f io 
Lr <_ 9 
LIVI 2 FINI 
Gen SOO WWW 
Ts et et fee (/) fe 


SD ad ee 2 a 
2 vtqoaww>O0 
AawsZoOo>rrt dey 
ZN SS wOOoeTwW 
ple: trae 


Gngew 


UL fete OR OL 
wote>sF oO Wea 


(Ske jet em 2M 


iw wItwWoOO 
mY <itrr CeO <_r 
SWE O eS 


— J) QO>5> 


no] w2JoO> 
2<I>~ONNO Tet 
et a NN = 


=> 
< 
= 


D 


¢ 


THE FREQUENCY OF OUTPUT DESTRE 


INPUT INTEGER CONTROLLING 


NPR 


WOU QO OU OQVOQQ OO OOD DOQQ OOO OO OO OD OQOQ0O ODO QOO0OVUVO 


206 





eleolelelolelololololelelolelelelolaleleloalelaeleleloelolelelolalolelolololalelolololelalelolelale) 
PUN DOPE BAO HNN FUE SRO DDS] ANMAPN OH DAO SAUMO FAVE ADO AAO PNOMOROe 
age fe SP PLO LS LAL SV AL AN SO OO 00 0 0 0 00 Be Be Re eR el Pe i OM ODM ODMDMMODOAOMO 
_. =\lelelelelelelealelealeleloelelslelelelelelelelelelelelelelelelelelelelelelealelelaleleleleala 
<>< KOK >< OK OK OK OK OK OK OK DOK OK DK 9K OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK 
meee 6 O50 044000000084004800000064446000040607000800404800440008,. 
>< >< >< OK o€ OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK DK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK Ow OK 


ODONONAOONAADODODOODOODODODDODDODOODODOAODOCODOVOAODOMOO 
Ut! 


oS) 


O 
T 
E 
F 
I 
0 
O) THE SAME INFORMATION 


R 

R 

L 

? 

N 
oGT. 


Cx <f Ow A. LU 
2 mt Oe 
ALAMWMaZDs 
><) pew et (me ore 
QO = 
me jie S <UL uy 
mm MCU O Lr >a 
Me YUU he 
Soo. ce j= 
Za > UwW2Z2ato 
OA CfoOn—2 
aqiu>a => Ou 
m QZ PaO a 
ODx+ Wish 
ONTO Oemw 
a OS ODO 
CerrteDaZzaae 
La wmmaZtot=S 


a] ~ 
TW a 
roe iL = 
(a <4 Lu 
op = Qu 
eS wc 
LW Li LL ae) 
Esto AZM 
<a) zac 
CNe Ler 
Ww © W) 
Ze2z Oa 
Ww<e W>D> 
 ¢<¢ TF = 
Ww Ow OZ= 
wd) Ti 
Om >< UWitlire~ 
—~ OW woz 
>< LL 
Tw xyYy MNwY)W 
Ld tel 
LOW => ng S| 
as — © 
= Ore<t <I _J 
ame | > Yan ho (7) 
=z ZL ae <t 
wo «ole 
S) paw Wc 
mlUCtiU NA 
=a D> wD 
poe Qe 2S med «6 <L LL Re 
Zz — 
=O > MWC 
QOocdw jp > 
a UUs © m=O 
eS She Ze oS 


riot Oh & 
Pm LL OO LL TE ert) OY) 
— ms | SOW 
ZFlmMaca Mum << 
— Ter Ol SLL 
mets (5 Te 
Wy ILUWSZ lke © 
mee, 2c CO 
—U<f hk he OZ 
home r wee> 
LLihiiLOu rr 3% © 
Pet eT pemey Net ee I LL 


Sik ie ai Ge 


HNAP WO 


SOLUTION, 


D. 


= G58) 


INPUTTING THE UPPER 


A REAL ARRAY DIMENSIONED AT LEAST NV 


fosfosleslesles)eSleslasles| 
CQ ow 
LJ © lees ad 
ae Fa) Olu —_ C) 
q) UJ jm § =e CT) UL 
QO tt !- == DN ad ZZ D> ee aw 
Ww<etl WY (a i Sb) YO TM pw 
=> CkKw Lj my ms) OC SZ LL =) 
Oe ND Ohm > (i ee fa Agee eee 6 By a 
WdoOwWW =voa WwW <<(iuww met LLY 
—_ ™OLU = tj Ww <— of 
<< ae z2a>D> Tr stu 
LL fame LLS fone = Ot> ell) 
Y J TT =>@ CSUs = Ly Gates 
wa *f CO Oz om) Oo =z ® 
ct 2@ CS % WwW > em LL) nt oe 
me <t JT <t ZNM fe JO WwW =I 
i << i aot Ve =e te 
— ZU) et ow Si =. de <ft ettY) O <i WJ 
— Ea Ore > he © <f wx _ 
Lnxrt = iad Zumt je SMW 
Oe = = MNxXZ ZwSseZz - Owe 
=> <« CL LL mre EMOW My 
eal mee am mt J MNO Wow ann 
CL” UL CD tet LL. Qk Iwas J {I< 
Qa==O = Ww Owe Yew 
= S=ew wo ad —Ot Wve 
morn, Oo) et) met oe Se CO TL 
ZW J — Ke be 9 oY wpe 
YW eH eo aHO «© See Wut te © 
Wwr2OAZ>nM Wee ec =—CO> Zee 
rT Wes MSM Ze © QOOMY bh °* me a VY) 
fm lL) J mt SPY tly Uc > Lu 
THIYODZ Ween BVWwWWewiLor JO. 
Ze <I, LUZ ee OyZwoataat <azZzr 
O NWe Omen ee) Ze ZU 
iL Ore Wea > WwW J OA 
jos LUI th em ee OLN 
m— OFflD—= wT O Wet Zee Soe RON 
= OKO _ <= Lu <I Az 
—— «608 OQ we mC tlhe eZ St Se 
aZnow ™.o et OO W>tnm Lu 
CO wmoe>eZ Amati) <A ooo cw> 
CS | | Onl a A ae OD Be 
Qe ON AQ ey SZ MO ee Let 
qg Zot OA. Cael OOZD => 
re JY Rt oe Oe RL 
= DPD DM Ww *Mw <hh a ae | ai) ol 
OOwdeNO o> _- 20.8 a Ww 
ZITWON Swe YO CrmUMN 2ZWNo 
mt > SO DQ <Q) ere wu —iUw 
(eS) eye a ® ty Sl oe = he 
Oe me WILL Ze) eRe O Vet KYL 
Wh JO ae 2D mm TO wW<re 
OTD ame Ore ee Lcthke) > 
Lip atee <tc WO BF DEzatDe Ww 
— Wazhr ol QO Je COW Ub 
2ZWwWO eke NO Vodadwozrt 2r™— 
memes DOYLE ED MO SOE ES 
<— 
_— VY a. 
z (a x< _ 


XU 


QQ OY OVO YOO YOY QOYO OOOO OO OU OOO QDOQOOUO0 VQQYVO OOO OOO 


Z.0r 





DOOCO OCOOCO SDS SIBWO OOO OS OOS SQV OOO OVDQVOOVD OO ODDO OVOV OO COO0000 
NE SPUN OP ODD HAI FUND MD BNO tN OFLU OP ODS] HAA PN OM ONO SAUNA PLO MOO 
DDARARAADO OOO OOO OO O ow rt rt rat at at et ot at tH IOI NI OI NI IQUE FAFA CA MOA OTE 
OD DOO OOO Oo ees ered eet rt rk erred rt md ed rend eet and cme emt med cmd med cmd ceed mm com creed med rend med mm ced rent mel owed ed eed nd emda end gent mm 
><S>< dE 9K OK OK OK OK OK OK OK OK OK OK OK OK OK OK IK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK D> OK OX OC OK OC 
nan000400044000464006,40440040000000000000000000400000000480 
>¢ >< Oe OE OK OS OK OK OK OK OK OK OK OK OK OK OK OK OK OKC OK OK OC OK OK OK OS OK OK OK OK OK OC OKC OK OK OK OC OK OK OK OK OK OK OK OK OKC OK 
AMDDADDODAD DADA OAD OODAANODADADOAAGDODADAMOOAODaADOMDOoOaDaaaOODeg 


” 
= Us ' 
Re Wt CW Oe CO ue 
aoe Ss Z7WOS zn me 9 OO us 
Ow OO tet oS Com LULU Lu uJ 
OQ wwOo@d™= LL. (aa Czre Li cr NM 
>< —<z7” = me © SW bm | wwe 
eMO We zolWw Use Z ww SZ — (oa) 
—O LmO<c>D > Us Usus << Zen) 
cS) ben U9 SS J mt =O om UL. pont > = 
ag Spe aor ¢ <_ ~~ =~ ee o<{t eS) [ai aes | 
JO Ou Oa Fo We NQ «OO zaeZzo 
Vara Wwe QZ wow LJ LL >= =>Zo<a 
<~O —x< <fh =S>> Aas ODOktZ wMOaZ CC LJ LL, 
Wits etl = oo D.= -IKe Cw Te Oo 
eo == NTA ‘on ) O NeaD> Www SZ bm <f 
ma SOW me OV eWeek SoZ wa> 
Ore AcyjO— Wwe Wi Ws WwW >< me Zam 
a oe <1 cm KO he wo MH bel Oud << Se 
Us<g merle RKO tle Zz en ret OL OF 
— ~~ e oo <5} Cle JO KL — —  <f 
~-~ > tnZ Le Ow Ow <worz (_) pond /) pons 
ews 2 e<WIO 2 Wwe NSCS SZ ig ae | 
LL) jew uy =— O Ce «K—*wOtY¢ li wD >OO = 
at -§ JON —r LL Ww WOW - LL SOx 
ao NaoZza> Pre Ri eRe Sol ol — > 
Ime IeCt>i %t~> 20 etvo Ch e& ly me ef 
pant med LS bee Jat = WW <f Wil a => LW 
LN wake De\M <i OB | Ww te OTUICOWW mr DY) 
< <Ix<OQ Wm wWZ Ze OO J Dd pao <I ns <f 
->W -F> ZVOe nO wWwwi<t wt UO = 
we << DL Ws woltl Aw LMNNO <a» Wid. 
bn <f fam be oie J rT OW Of Usa Zwudc a Gse ae Ud 
= OZOZOQ We eo wesc WY CS CO LL tat pee 
wn WWwoOowWWw> > YM DW"MwW—_est Se ee 
all aan I = a | J <r — Tale) —~O 
235 OZa=s«<O <ke ey WOWTAdA Se ri ae: <a 
Woy OCS LU > Zam uw DRKeOQOd< crear = 
QOQq MAW > Oo 2s bw J dott) © —_ © > NJ 
YW> Zw WWM Wel D> att wwtnm mW ~ 
Q waa cr Lea ee Zt oR oor ww ONx<Z «<M ad 
Ze S22 O RKO We OW XN LsS~w <fiij~w 8 
mnt et eg LL LW = C™ WLMNSZaHQOOZ << ml ay LU oe 
a © weWEe nO Co + OaW mL e<tiy — 
ma TOZOw -—aQ Fm LL) mt bee ee = = N 
Cm we ZO ed nm Oe ODKRZOOH IW —ZNcCae 
< qq SW Rw CK ZVeoOr azz ZOenWO OD 
Lu MY mas aY C2 2 202ZZ=>0 mat mt Se et 
. w os act WA teDicetm= BD =— TK e 
Pm QP mevy eS OO KX Otbauomoud ct WwW Ws) = 
Oj Mm Oy DH ¢<o oe Oo wz Zz + 
> 3 FP w © Lt co Lu f jD5«esS lw 
O—= foO—WwT e 3 ~— — MM QlLewtpiL, WY 
Zaq wa7ae>rarowo NZ WEF | D0 amNMNO > = ww <x WN 
2Y~ CDC WO™= =O feO it Tf} Te | ad Zt Jj < << NOW O20 OO 
Oh Peers TC Ze iL «<_ < ty e@ eNeeYese eS 
On fora wt mil 0 W2%N%527Z 2 DaOrVoORK-O0O MoO 
vi ee © me - nT ap) 
Ly LS ee Le RNS usen AU 
ae faa a  =eOMNeH M pet _ e 
z} = Li Q Low mK De men ee II 
< > — = welweeoe UO N IN IN 
< = ewww e~~ = t 
>< ~~ LCNNWDaDdwW ys paw 
us O <t<«Ka<x<xS=x<xK eZ 


jawe o.. © 
Y) =) —J 


OW WO i) a) 
QOUYOUQOQUY OU OY QO UO YOU OU OO OV OOO DO OO VO D0 QO OVOUVQOUUO 


208 





aa oe) OOO OOO COO OO OOO OOOO OOOO OOOOOOoOoOOoOOoOoOOoOOoOOoOOoOOoOOoOOoOOoOOoOSo 
D mt AD POPP DDS MA QO FUDD PD DPS MHA PUD OPM ODO HNO FIA OM OHO SMM SLA Oh 
FPS ne PEST EOL LVN ETT LE NO OO WO OH WO WO ODO Pm ee ee Be Pee he OD DMOOOO DO 
ed ced ed ceed geen set eh pend eed pnd ge ed eed od cee) rel oe) ol ed = el eet md ed) od J Sed ed oe eS eS et ie I Oe es Se a Ee ES SS ee 
>< >< >< >< > > > > > > > Oo > OC Oo >’ > o>’ > >’ > >> >t > Ot >>’ > > > > >> >= > O>*’ = > OO >=’ >= >= >=’ 
eaeao. O.6.0.0.0.0 0600060060046 04060604006046006.0446400464406000040464409044 
<>< >< > > > > >’ OE > O>= > >’ D> > > > >’ >’ > > > >’ > >’ >> > > >’ > >t’ > oe’ oe’ >’ >’ >’ >’ Oo 
DAMODMDANDMNDMMONMMDADMAMDMMAMMMMAADMDADMMMOAAMMADMMAMMMOMMMMMOMDMADADAMMOMOOD 


~ Lis 
-_ a 
< ’ _— = 
uy = em LU 
ad e NZ mrvach 
° MmOQ 2 | 
on) wood LL Us er lL Wn 
=! bee I= 4 hee Nand 
LJ pm <f jo ON SL 
0 me <— ad <— CL ot 
PiU mm CD Le th ee ee 
0 rm tO OrtOr wm 
—_— 2 ow “L=> re 2nd 
cy o a aS HJ. «© 
ij <= AS) Ors <O pe OC 
ome =i ah me = LUG ef til 
o~ CS Out a 
Za mie OQ oi gel et b> ne 4 
= ~wW = fa TS OOMOW 
~ N= jm mae ~ ORrOOc > 
= x pang CE al _— ZeoZzoOa 
—<~_ we Woe ioe) <2 ftow— 
x< & ant foe es} aa] Or CoVWOR 
oo ph JY = uy ZNn<t ze 
=> Te Qo Ne} <_ HZ — 
< ‘aa > a) i wa Om 
er Oe wad ae) e xm eal OS 
: QO Wu Ye) an Tao ban TOO e, ee) 
— fp o& Or a) a a = O00 UW 
on <f o ze e ~ ST me SLI OL CL) and 
maw =) We) ~ 2= aI 
<WOW oe ee WS SUS Oa 
oe ee pond ad e aN | m=<LONAN Se 
cl A. — + >< 08 — 
eZ NW oO= 9 ott LID 1) de ALU 
tz —> i} — ® aon OZ O— Tt 
D> ted Lis Oe Zeus The 
Za Met ~ > — mone OY Iazar~e om 
ue 2 > ral CSine- a, — | meet) = 
CL SP eet jan << . —~ wo-= MNCLO mee Zz 
a “OZ Lis Nx 2D ae m<—<t iP OOoOr™ 
Zao WILL. < Nhe > ~e UW 2 Ue Lo 
~ it = _ % ety LL, >< Se Cl mt LL 
>— LJ he “iY i of -— Oe Ze FO 
ZX eTWO) om LVN® LJ ¢ Tij—Y) e 
aaa YW Il CO eum > e _ 22z22Zzwn”N 
e-_ > fm Li NI wu meg ON Lis Lh. ee CC LL te 
> ns Zw A OD ed FOZORenL 
ay —~ ~ONO —_ €) 4 —= (SO 
“~™ > oe LU! >< <I jaw oe e © > ee xfAOMw xr 
~~ & Ll wat OF LL wf eat eft @ wp “MO LO OL amo «3 
~<a fe LL) ee /) oer bh oll) Line Ode pas e 
—~a~ ™“S ~n x + jt iLx< Or OL LU LL de Y) 
om WO il @ ee, ¢ —~ = Ud WL A A 
<a—N 2 OOZa—— | = Zz Wan OZOUDZ2 
LAGI m.0 Crt ee wee mom > O=— ht xX< <M 
NO me |e ma UOMNXK KO Iw He mm) ~a CO wwi¢gerze 
it lL <f fm fw LLJ << oc mm WI QO TeWwuw<t 
un ate Se a. C)<fpe LU I) Wi) emer omer = jj Qu — BD LC rea uakLawWN 
Iayve OO a = Os my + bow Pr i—-o 
C£2I>a YOLON -Z D wt NL ere LULU D—= Nuss w 
2Z2a— SuL~N MU LL LOOK Se OL UO mit Ol 
ax > eles 
=~ LLJ=> vom ~~ 


O O OO OO O OW OO 
VYVVVIOVOO OO VQOOO QO OOO QO QO QO OO 0D0 OOO OVO OVO VO OV VOVOVO © 


209 






Io) OODOD OC OCOC OOOO O OO SOQ OOO o COO Oo oOoOOOo OOOO OOCOoOOoOOoOOoOOoOOoOO 
ODO KAI FUADOPE DNS AAU FIDOQP DHS HSA HPO HD NO HeqaumfNVWMHONS Sanum FN 
DOAANRAADAAAANOO OQO C000 OO tered cat el at at et ot ed EOI I IOI CUCINA MM 
queen pemet peed pened wend nd geen emed peet coed weet OI AI IO CU OIC NIC OUQU NI QUQU CAI CU CQUCQUI AI NI QUNI AS CIN QUIN AIQUAQIQU NAINA Al 
>J SOK OS Oo OS Ov OM OK Oo, OS OE OK OK OK OC OS Oo ORO OO OE OE OE oe, OKO oe’ oe’ OE oe’ OOo’ >’ >’ Oo Oo’ ore’ OO Oo OOOO’ 
6.4.6.46.6..6.9.900406.00060060664.46,46.4.0460.46064000460.40064464400060.40440aaa0, 
96 >< dS HE oe KE oe OE Oo HC Oo. Oo Oo OS Oo OKO OOo OOo OE Oo OO >, OOo Oe Oo OE Oo >=? Oo oe’ > Oo oe’ Oe oOo’ > OO 
ee ee OO OOS SSS SSCS AV ASS SSasnaaoaon 


a) 


LL 

Qa Nn - 

wn 'ae wl 
oe ee eo ed * 4 
mt) IZ a<i<f 
WWW LL <f Of wat 
mat /F CO wJWY) 
*-O WWM <a aw 
Ym NAQm mL 
my ODCMWwWnom 
Za Cea TE Oe 
— JO Hwa 
IM) bh Wee WWW 
CWO eZ Se 
cond eek on En) en @ 
Nat mi COzZ 
2ZW>= QUO a 
© WIL LW JO. me Lh + 
WO COWJO= tae 
YO. tOMN@ Fe LI He LL 
ce I Le EE Looe DY oo 
mati £2 OO) LL Zz 
QreOQN et Twihataiin 
pod mt et LL fem ST dd te 
wJ3ZOOWwW &- O 

QL met baee oe Fa = 
= DNWM ext oO 
mee OQ ee Se he LL 
C—eTOZz KOO! 
et Fh LUI Oa lL <t 
Sa WEewyszs 
RKeO>raWwOoreO © 
Zz Y~Y Qonanm 
Sww<tN fla Ze 
srt NE wi 
PreK NTANndtox 
< CLS pm LL <OLLS fee 
a eZ WL f& 
COOWWOx” NS eae 
mel ye ie LL Oe LL 
>ereuw COrd S24 
<I wiewetiio WwW 
LS ap LL Se te ty LL A 
So.) Z2Z>O ete 
Om wt Woe 
— >< (TNC Ywes 
YD KFA QU>FQ 
YU Kefl LL 
COE wes CD gy LS Ly ST LL 
Agro ee Oct >> 

aol aed? o> Santon] a) 

Cus NOS Tia bee » 
(on ko 4 ee ee) eo 
LO=ma> iL ZO 
= Ow tO Sw 
CO wawQa =nwSa 
Weadqea Oc<kKO 
YO meet eyo 
OO thew rom <2 
Us joe Oo tL SL ey 
GZae2Zworwr7aa 
QO mw OE OL ee << 


S} 

a LL 

YM SO 

Ze oe 
OOO rea 
et Po 


pee mt LLL ey ome 
aikOr->Zt Te 
meZZe CO ekN 
WO Wi ce be 
aAaM~A >WwoO SO 
=2OnN FDU 
OD Wee >= ae 
QOUOeKTO ode 
Ct tw 
OworewJYT ecZz> 
et a eee oe 
Ne Ut COWL 
I —mwH wad LL) 
DBOWEtLOaAZOT 
ZZ Ss CO LS et CD fe 
Weeoao, © 
TLYwowwo 
m Wen OeZzy 
Nas Ome <t iy 
eZ OAOEMN <_t 
OOw ake VO e 
Wes Tet zrOoamnw 
Lhe Wa /eO 
NO PO O- 
=I rZz MN CS 
ad OL OLY me OL 
Drewes ms OO 
GINSseRO Kx 
mz oO Were 
Ne ox ZrOO aa 
Us ts CD K<o e LE 
UL fer LL Sad 
Mw OaAUOW 
ret LL LULL Or LL LL 
T>>rltUeE2s 
<——_ oh 6hLL COS 
Ws We ULL 
aw DOZAcax 
Oo eWwC) ! LL 
=O —Ovr te 
Ow Ome Wwe © 
QMO >ouw 
Ie WK OZ> 
awJemme tit | LL 
fee Tihewn JO 
rt ped Te LL CO bee 2 LL 
eZ Sle ZO 
mt ep ee Of te 
z wile ~ oJ 
me? eth WW Oa, 
m= >> TOW 
Z2 WwOore<ngw 
Zak oO. << 
Od Jit. WiLL 
Www <fte i eT 
ret a > > et eer SF CD 


=J >< { 
OJ LW 1s 
ek ee Ee | root 
aL eS 
SMNWWiNeWIo<T ie 
“OM >OQODIwetyY= CO 
fame foe LLJ OBI IneaZw 


M OOWQ pre = 2 
Ye QOZTORNYO Om 
Ogawa] BVWwSZF ma 
ZSFENOQ NWO °O- 
LeZearZze 2502727 = 
NID>WOoOOCat>arKorat 
Zit Wwa=D> rot tnt OC 
SsDWLORK Oto Seat 
<q O0@0QO OW OFfD 
=Oae DwWweyw27s 
CLL Sot ee Te COLL 
Sd LD ka SLLeem 
TWeww ZVunz — 
YCitIe~ Ct ame We = 
Caddo t-$Q 2>M 
YM  ZOCIO OW et 
Lime Amel i= ie 
SHOR RU >NWeOmW & 
BO AOkeEOW We WwW 
SAIS MOvw Lz 
Qt IJZO ecto Oe 
SE dh Se Pe ed 
wokwWw a<aAW<tnt we 
water JS We oS 
<{ MAM LYM w 
—- Dw Je O= 
Ow ote Z2nwy fo Wi 
ZID CWT emLeE> 
ON? WI LL Sot NO OO 
Nea O>D>WSTKENs— Oe 
ma ICE LC mee lL  e 
NNSKatwaZ2o0wsrt= 
Zea NWA =O> UZ 
OSU VILLI OL er te 
—“y ~—~oOn Zoake-> - 
pmo LL} CY —JI<I 0 NZ 
OaAawWreeRKOAD> wOowle<f 
WtreOFwm WwwnireWr © 
wd bm SLI OY LL fee Ce 
LON JF Bi 
ws mt eA Oa ee 
fad baw ry? OD LL Sr) et fen 
IO *fa*#t#oxr~O eae) 
cus AD rfoOw <= 
Ly Or me Tt eet OY 
>WOF Zw Ti op) 
DOO vLeaS<INnNT>o>atwed 
bond frame nd rt pws TO 
ZwZLxr = eo wW 
—\7w OO eMCLO ~<_) 
OUW*SQ0QOwWwWO —_YD 
ZW CO Te 2Z7wo 
Ww Wey Om ey 2M O 
TGTNMTS Teas wSaTOde 
TSS at pe LL fee per “SS OY Se 


REMARKS 


QQQYUOQ QO QUOD OOO QDOQO0 OO OO DOU QOD OOO QQO OOD OOD O OO QVO0 OOU OO 


ANNO, 





DOODWDIDOWD DDO OCOD SOD DOO OOO DD OOD ODVO OV ODO ODO OOWDOCOOOCO009O 
DPR DDS BHAI FID OE DD SD eg NO PLANO BO ONS BSH CQU MM FIDO DONO SAIM SIN OPM DPO aANUMG 
PAVED CO) SSN Se it SPLAT LV LAST LSU 9 0 0 OO 0 0 WO OB Re fe fe fe Pe fh Pe fe OO OO 
CIAIQIQICUCIAQIAIQUNIQUINIAI QU CICQURQUQUANY AICI AIG AIRING AUC AI CICU CUNICICQU GUAT AY QURAICQUQI QU CUCU CU GUC 
>< >< > OK OC O*K OK > oe O*K OK OK OK OK OKO OK oO KE Oo OK OK > OK OK OKO EEO oe” > Oo OK = > O>*K o> OK OK Oo OK Oo” O>*K OO’ 
fee Oo. 6.0. 6, 6.6.6.0.66..6.6.6.4.46,46.6..66.6.6.6.046.4.666.6.46446.64.464.6.464460464406 62, 
dS > dS OS OK OC KE OK OK oe, OE OK OK OK OK OK OK Oo OE OS OK Oe OK OK OK OK OK OK oe OK O>*K KE Oo OK OK OK oe KE Oo ot KE OK Oo>*K OK = Oo OK’ 
ODMMODNDAMNMDOMA AD HOHANDADMNAD ADA ONMNMONMDNADDOMNADNDDADOMNDNMNAAOMNWOOM 


OO 8 ect 
At —>Tt 
mt ot J) 

— CF Sie 
WiLL e<fOus 
CZ Gao 


fam ae  LjL 
Je OS om pe DL 
wae fe ot LL > > “SS 


(<fae O<TEON LO 
> fe LL ZW J 

ZO. Ub Jct a 
~OFZALNODwW 
CY wet LS 


me) CO OC rm ot fmm fe pee 
Oe 2 WZ 
CnwgTti WwW 
WAOoZFo 7a +O 
Fr CZIZMN 
fer TO UWI Wl 
mm OF; Qi 
fae me Se LL = LL CO — 
IFK Qo aso 
= am ZF Oe 
ae ee ee a 
i aa @ rome Ge 4 
Cw7>>-> oot 
wajCcasoaa > 
— WW >» 
00 fhe uUee 
Ze oe he JOS 
2 2 —Ww 
WOd => ke YS 
Mis  OZiuism— >] 


a ee aa 
me) mt SILO 


SUBROUTINES AND FUNCTIONS REQUERED 


ARE INTEGRAL PARTS OF 


‘EF BVS 


AND FUNCTION 


LL = 
~ 
<x eO - 
— om NZ 
Lee = Ow 
Mette = 
a >< Ld LU 
eDO NZ « 
m=O Met 
Mu mW 
(a Mi ae age ee A ee 
— i fawfee OK 
Lhewo<t NOOO 


@ f= LL 
CL Le I oH LS 
LL) et fm 7) a me OO 
NE CC eluvw~> 
= CZ LL) me << 

fooe fore (1 me UL OF LL 
LWLIWI<T At fee 
TZ 7 UL 
—H-OJDTK NM” 


Othe © Oat pee 
~ > ® 
OU) wl *t_J 

= WwsOOO 
Q04Z2>+ €£& 
LL = LL C= 
m= JED euzZz 
aA <P OUOOO 


6. z= ~ *1i) 
CO. LL 
=> 2dswil) 
NWOAOnNna 
Lees Zz Oct 
LL) fb <t <f{ ja 
Oo OxYCMDS 
WZ DGe 
fom fom “) LJ 2 a= 
MAtiiuwt 
2D Le "WoO 
Sed OY 
qT 2 Pu 
S) > fw SF tt wee ST ee 


= LL CO OL 

LL ee 
eT Ed Of) 

ODO O8ENZEaO 

ZPnWe=Vr200Z 

pow ee OY) SO St 


REFERENCES 


METHOD OF CONSTRAINED OPTIMIZATION AND A 


"A NEW 


Mao Jes 


BOX, 


QUYVOY QO OOOO OOO QV OU OOO OOO QOD VO OO OD QQOOOY VO VOD OOO U0 


PALA 








DODD DOV OD OO DOI VD OCOD ODVODAD OO GO OOD GVOQDOTOVGM OVO DOC OOOO OCOOCOoO0000°0 
FU QP COD © St NS SULA OPO BO MAMA SLD OM OAS ANOMALY OM DAO aN MPN OMOADS — 
MWD ODMDMODAARARARAMN HOO OOO OO DOO wh dh et rd et meat cant end wet mt QI AIO CIO OI OI OI CO OFF 
AIQICIRQI CI RQIAIRQUICQIQI AICI NI AIAN QIAO SO AMA O MMO MAM MS MM HAMMAM MAMA MAMA MM MMM MH on 
>< 3 DS DOK Oe DK OK OK OK OK OK OW Ov“ OK OK” OK OK o>’ OC >’ OC’ OC OK o#eK_ OK OK OK OK’ Oo’ OK D>’ >’ >’ OK >’ >’ Oo OK Oo Oo OK Oo’! OK >’ OK’ OC O>*’ 
O460600A4440444 4064464404404 4608400044 46644, 0664 424400 466040604640 060060 60006464 
>< oS DK OK OK OK KE OK OW” OK OK OK Oo OK DOK OK Oo’ OK OC Oo Oo>*K’_ Oo’ OC’ >’ Oo’ OK’ OC Oo! Oo! OK OK OK O*KE OE o>’ OC > OM Oo OK OK’ OC’ OK OC’ OK’ OK OK 
DODDDAAMDADADADADDADDADAADMDADDADADDAODAOADNADA DADADADAODDOOGcaAad 


6 a 
e ee | 
= e co 
ut e 
O e o~ 
& ® CS 
A e > 
e =a © = 
ce << © ~ 
a. 6 -~-~ DD 
< > ad « (oS ee 
‘a al @ LJ 
(oe) a | s) at o ‘am 
Lui ® - — Ww 
~ a. © =a > cS 
Ee) — e Se oe Pad 
<< 2 ~> = < 
_ ee uy 8 ~ VY) ae 
fa Pal ~ e@ 35 3¢ i) 
a os) Oo © co za 
© = = e ~-  & aay 
= oad ~ ° => = 
< N © ao w = 
mw md © ewe A S) 
Lu a e oA — -) 
- = > o = 2 = 
=~ (a e e UW = —Ja Oo 
a. _ < @ WwW WO — nn <t 2 
=m C.. Ud e >» C) aad _ 
© (am — ° ~ ~ — <2 © 
>) = =) ® StS = — tu 5 
—s e cm W! oO. a Ut 
~ ~ ia ae @ ~ AW = — co ke 
o e ® ® i hand fad 
WY we rh e —_- = Www <_ 
a} ) o™ e a e Y —_ << 
© io Oo 2 9 ~ YY) cy —_ (oe 
aa LUO. — <f ® (a ~ =) = 
— jouw oom ™ a ® ec: ~ e Q. Ww Uu lL 
LJ nag Lr © ® z=z_ = = YS) (sa 
= Wo a e -~ Oo it — ht — U/) 
ULI nad Se 6 e > wW fm al — = Ww > 
fa l= onal oO ° <i! () Cc = - — ea) 
Lui C= orn ° a ee Ny = — @ tls — 
am NX Oo > e - D> j= = a Oe Gh. 3 oO ‘o 
= t C= @ e => WW a = (ot nod bom a 
en) a= mttLLy OY ° z e < Se Tear 
zat ~N=O0OO ° _ * | on oo yo a2 2 2 = © S 
ais P= qa mats) CV LL ® om as) < VS) 5 So) 
_ Ge) =O e x OO ra d a - me e e 
—_ o{_) SANZ 2 wa OW — nm ww © oO WwW =< — 
= = we wl e a. o Za ® wm ZS BZ We _ jwe 
5 {mw te ® <x Oo ag << > > 
ae vad wo x«< ° a Wy -_ © > ~— = Se e ® 
Ow Ws aw lh. CUS ® ao - © => PL ee a UJ Lud 
Min Oud — Ws e > e et wg w WwW ww W> WwW cE 
=) CO ¢Y —=OArKGA e Ly OOm OF tw * ete Fane =~ 8 
CN ate te Tl e oe 1 OW o> wm YY CEC SC WL etm mt 
{fr wKH = MWY) ° = CC woe Ws OF D> DOD code et ee et 
G: my) = © ht CY” heme bd — = eNN JZ FP OOUW L& YES) meee td ONS ee 
=e >t ft Yoer> a iy emt been © (= It moe 31 De 
OQ we fe eWOW ° (Gls | ies Ae | lt oa) Oca 
Oo aa OF wrwrudw ® ce WwW Wot — it Y Plime dt it~ It 
© ® oOo =m <I Wot jer LL 
io e aS am > 6 (OOoketL COL —_- Oy Ff Cm UL, mt ee CD LL 
C.. r NM QZ SWZ Cue ze eae ae, CO ee ee > CD et tat 
om) 
ond 
OWOQOU QOOQQQOQQOIOOU0 OO =) CoCo 9 “Wa 


212 





SSO OOOO OC OC COCO OOOO OOOO OCOD OOCOOCOCO OO COOOOQOCOOO 00000 oO 
OID ELD OPO DO HSNO FIA DOM DDO SNAP IN OM DODO HAM FIN OMDB SOM FINO MOD 
rataaTadtad Laalaa lad Loa he ake abe h< cc ec hs wee ew gd oda atlalatratliallatallaopeloneh ener enelenrb int. ini. ie tet. 
AA OAD MAO OO AAG A AM AMON OMAN AM AMAA MH A MAMN MAGA MM OAM OAH MO MAMA MMA MAC) 
>< >< de >< OK DE DC OK OK OK OK OK OK’ OK OK OK OE OK OK OE OK OK OK OK OK OK OK OK OK OK OK OE OK OK OK O* OK’ O*K OS OK’ OS OK’ OK’ OK OK’ OK OK OK 
666464660460406046646464646464606464.0.46.4644446,4044644.40406006460440460460424, 
>< > > OK OK OK OC OK OK DOK OK OK OK OK OE OK DK OK OK DOK OK OK OK OK OK OK OE OK OK OK OK OK OK DO’ OK’ OK OK OK OK’ OK OK DO’ OK OK OK OK OS 
NDAMDMOADMONOMMOOMOMAAADMNMDONMODADDADODDDADDDDMDGOHODDDDDAMDOADHDODOOMH 


pt pa | 

feels! 

eo Cae 

MW) 

— ma 

<< 

~ % 3 

ey) a. oe. 

~~ LU Ly 

~ ~ & 

0 ac. 

wat LU LL 

FP wee we 

Ww amt pa 

= Ox < 

= pmo ef <—f 

we == 

= Oat 

O+- 1 

— Sy 

© se as me 
Sf al ee? ae pe 
apap > eta 

mm Ina 
~w Ot} if LU t 

ie4 eu i 
OO ——=Q,, — 
Ze meg eae pent cee ty prey 
[fet Co meet peed end nee 
rm LL, we LULL J > = 
>= > a 
NMA i 


NUMBER OF CONSTRAINT EVALUATIONS 


I 


NCE 


e 
LL 
— 
< 
Zz 
(om 
= 
fae a 
LL 
— 
Oo 
VY) a | 
Lu -— 
‘ail 
fo 8) UJ 
< i 
Comb 
oe a) 
< LW 
> cS 
Pad 
LL < 
ce a 
WC) 
e = 
Q = 
=z 
2 
VY) ® = 
Li}. CQ — 
= LW = 
— x < 
=< SB) WY) 
= es | 
“ ee | <{ 
<— — 
tt fa a 
VY) — 
= -_ 
aa = Lu 
a LU => 
= = 
WY) LL = 
Us > =) e 
ey <IO WO —_ 
— wie LL) + 
_— a+r “NZ <—t 
(a N> 2+em eT 
Lu =~ Os ec =o 
> C13 ett | 
> i. v<at 
in, ot = fl LL. 
ox On O tt 
NS i Sn || 
~ C= et Th << 
LU a Loe iO OS 
om Od MQOAYYKW 
SY SZ SSALwLw®s 
red =) =) 
=z Pa Pith 


213 


) 
&) 
Ze 
a | 
| 
ee 
e Lu 
or) foe! 
(ama) 
‘os -_ 
CO 
Y Pa 
om. VY = 
LL ~ 
we W) C) 
op | —_ 
— —_— (ama 
< (ed _ 
CC ms WW | 
wo > re 
Zz ww -_ 
Wy + wl ie 
Oo at <_< LLJ 
JS = O 
(= A 
LU it o— Lu 
oO 2m 2 eS | 
= Zw = fae) 
=e tore — 
Pah a. Tp) 
~-~ <{ 
= © ~— Lu 
om) o —— UL. 
CO @ tien 
=a WwW NWe LL. 
<t e oul me ad 
wm —_ et == 
GJ => ee? cf a 
LL f= = a ee *~ a. 
OWgY Us 2 b— ee 
Ler j Le 8 (LL ~~ Ss jf 
CS) eat tomes nee Aja et 
LL) 3 ae #2 tl tron a 
Udo © LL, © mod tod = 
YM —_— AS SS | 
= =— il and a | 
wy = we iO “tke 
ee er eae ee am 
> OW poe > DB OP Doel oe Lone) 
S) wns ome ire LikL Oba oC 
me Pas 
— Ne) PM LU 
OW OO >) > 7 | 





QOWVD OS O QDQOO OD OO GCOQO OOOO OA OO QGOQOO OO GOeCoooO oOo ooOoO00oo 
DAN PIN OP DAD HAI FAD PD HPO MAIO FIDO DNS SHAMS UOIOMH DDO HAC SIN Oh 
DOD DODD D DDDAAADAAXHAKGAOOOOO9O90 O00 we et at at ot ot SH SII ON OIAIOIOICUNIAI 
Maton Loa TaaTaaTnaT aa laalaa loa taaTaalaaTaataa ani aa] an anlage gk gl cS i Sk 
>< >e >< De >< > OK > OK OK > OS OK OK OK OOK OC OK Oo OC OK OK >’ OK OK OK D> > > > KOK OK >= OK > > > OK >= >= = OO >’ OO 
eee oO. 6. C.0.06.0.6.0.0.6.46.6.6.06.0,.46.0.46.6.46.46.6.0666.06440646.4.4464604460440049 444 
>< >< >< >< >< >< DO ove DE OK OK Oe OE OK OK >’ D> OK OOK DK’ OK OK OK > OK OC OK OK > OK >’ OK >’ OK OOK OE > OO > OOK > OX >= 
DOODMODAMOGDANNDMOODODDADAOONDODDAMDNMMAMAMMHMNDADHMADONDMONMMOOOMOM 


—_ 
- ~ 
~ un 
1 ? 
oO om + 
he) Ay eA, 
—- oO => bmn 
= wT ~~ 
ao =. 
Fae a CO wee 
nd << ww | > 
= ‘nl AJ (RB 
ee “~—_ 
© c% wz 
— er j ~— 
eC = AN ef 
=) CO t~ iI 
©) -_ it Use <= 
<< CO >< at 
= ere Cf aeJ eo 
= wu GF ots 
=_ ZOAm— OO + ww 
a) Witt L\ me > 
2 OuMmnm On 
© Ta) SP aw ome 
uy => @O0e eta 
O 2 Ue +m ° 
° = OveIYex< CW 
— st 2S egg Wty 
= Fo Daman > 
——— pO | COU mA. = 
= | Moons | I | Men Leeetomnd 
~~ OO oa wee joo 
Sect << a 
Le CO AGLIOK-(LO 
— CQ Samy YS 
<{ 
x foe) 
— WO 


I1))-EQ.0) GO TO 13 


NPT, NFE NCE yNVeNVT2VeI»FUN,CEN, I) 


Ge «= 


214 


10 CONTINUE 
I 
L 
R 


FEASIBLE CENTROID FOR STARTING. 
TO 7 


eNE.O) GO 


(VCl,1)) 


PT »NFEsNCE NV INVT »VaKyFUN,CEN,O) 


NG OF INITIAL COMPLEX. 
FIND THE WORST VERTEX, 


QO TO 1? 


‘J°TH. 


THE 


l 


1 


J 





DODO ODCVOCOOCOCOCO OC COCOOOCOOCOOCOCOSCOOOOoOCOOOOCOCOCOCOCOOCOOOOCSO 
DONO SIU SD OP DDO HA FID OP DODO eet NID FIN OM DHO Quem +F OOM DAO aim fury 
PNrovTaaran [paren Thales Kat lua las laa ts she mt hake hee et eed at al al alralallatlallallavekvolvs ojokverolelvelel. Tint int. 
PP SP Se Se te SS Se sr ee ee re et ee ee ere ee 
>< >< >< >< oe DK > > DK D> OK OK > OK > DK > OK > OK OK OK > > OK OK > OK DK OK > >’ DK OK >’ >’ >’ >’ D> D>’ OK >’ OK > OK >’ > >=’ 
Aoaaaaaaaaadagadcaaaqadaaaaqaadaaradaadaaadagaaaawtaaaaaa 
>< 6 >< >< >< >< >< >e DK DK OE DK > OK OK > OK Ok DK OK D> >>! > OK OK >’ OE >t OK ot DK > >’ >’ >’ D> oe >> >’ >’ D>’ OX’ > > D>’ 
DDODNMDAMODNMDADNDDODADARQAAMDMOMDAMMDDADDMHDHODOMNDAMNMMAAMMDAADAMAMAAam 


eGEeFUN(I)) GO TO 16 


K 


il 


2 


I 
N 


EQe-1) JN 


‘ama | 
—_) 
=) 
fo 
jue 
a 
LL 
=) 
LL 
is 
_ 
CA @ 
cw >< 
< LL 
= ‘— 
© a 
= Lu 
> 
=a 
< ~ 
— VY} 
cf fa a 
== ae 
= 
VW} 
(ee) — joe 
~ =i | 
~ Lu 
© —) 
- aa) UW. 
= LL 
© — ow 
co a 
=) [ea 
Ry CQ Li 
my ra) = 
CO = Wu © 
amd ees > 
=z > —J CA 
oa Fas O<jc Pl 
fm LL >> <f 
© 
Ow LZ Q 
9c) oe =—4 
e — am} 
Yvonne CE b= ee 
“= Ly = 
—_ om) mo = 
Ow W STDm— WwW 
mI SD Dw -) 
nA = Il 
Oemu, © Ww 
weer jf ee if a 
<= =) 
OuwuwzsO mie © 
Omm SD = 
mt CD © 
Tm OO 


OOUO WOO 


% 


HA*VT 
(VT+.5) 


= of 


Or DW tL 
O>"NO>— 


Zales 


ERVED. 


Sj 


CHECK FOR IMPLICIT CONSTRAINT VIOLATION. 


THE 


OVER-REFLECT THE OFFENDING VERTEX THROUGH 


TIME, 


=] » NV 


DO 21 I 








DODODVDODODDDOVDO DOD GWDDOD DOV DODO OCD ODOC DOOOQDOODOOOOCOOOC0O0O 
DP DAAO HAN MSFIA DOM DAO SNM FIN OM DAO HIN OHNO DDO ANUMALPINOHODOAaANM 
ml OO OOD OD DODAHDAKAAADARAHNS OO OOO090 DOO Hee ed cad et ad at th ed MOI 
BS EN ES ate IS BIOS RS de BS BI dO a ate DS a i NG ov Os eS eo ee a lrelralral allalialial alialeviallallalratlallaliallalralalra’ 
>< > >< OK OK OK OK Oo OK OK OK OK OK OK > OK OK OK OK > OK OK OK OK > > OK OK OK OK oe OK OK OK OE OK OK OK OK OK OK OOK OE’ OW’ OK’ OE OH’ 
0400004 00000444000400040640004040406000044004640604000000008 
>< >€ >< OK OK OK OK OK OK OK OK KE OK OK OK OK >“ OK OK OK Oo OK OK OK >’ OK OC’ o>*e’ OK oe’ >*e’ OW OK OOK Oo OK OKO >’ OOO’ >=’ OK Oo OO’ 
ADMNDAODDANMNMONMNMNNMNMNANMHMONDNAMAMNMMNMNDAMMMDANMMDODMONMMMOMOOOMMmM 


e 
cs 
po | 
= 
S| 
VS 
WIZ 
* Y/) & 
OW 
Q. Ww. 
CL 
a a | a 
e ah Q = 
‘an o © = 
— eS = Ca 
= het = — J 
[of © W) (eee 
— ce ~ z= MZ 
Pag — = a Oz 
uJ z « us a UC 
=) uJ Zz e Oo 
W uJ LJ OC <lo- 
(eos) ©) > UJ — =) 
fa Ce & i cS Ow 
<f CL OK z = Pad Pa 
= Tw = a) < “> 
© x 4" Li LJ as ue 
=a - (alos ~ = O LL» 
—_ b= LL ~ ae} za | 
~d -_ > ~ oe) = YJ 
—_ z o > © el 
cael — << e LU az ge 
co © tS} =— a Ww => Co 
* 2 => > -— Y 2 ole pe sf? 
_— Oa = I —-— We men 
DLO pent uu uy Lid > ~- uy <I OO 
~ Owe rm ® ~r z << uw oo bp sh 
mt he ~ CO ~ = ws NZ 
od anad 8 8) uu — thd — a C) ie Fe cO co 
>> > > <> oO <I Tm OULWNM A Li CD 
5: ed aed O ae Wat z > >Aa0 FSF amc (9 
<I > = Pe rk ICA ~ uJ waa oe <A Oo wn a= 
LZa— bap | Mie Ud Pi i P= = eO —- Om 
QC. at aed om LL Le > OV. Och te cali, 
Jaz ee met cf >OO 2 Oo ~~ DK See a) ee] 
<<. = =z ad az_ ~ Zz =~ FKFKLa =—OMNO Oo &— Ou—- 
ius Bs] > il LLJ ret“ pee So) 2 St 2 
— <f =~ = aIeO ef& =) ee 2aw Oct -~ © =Z 
=~: = mb Mt) CD b= Ls DeY jl eOm - (4 >< LL 
o> = < ont > rep its < » — tlie JI §-¢ WO ww uj ® 
— xe J a SS) jame ke <x > Zt CO MNZaZz —— 
wont oO  »>aa ILO—-zZ UW ~w (Lou time ¢ Oem CY aad 
>= — Zh ole LLL) © Aju few oon LL ent LL, ot Ne ete end LiLj 8 
te e<f > ~~) > LL. of LW uy OS oe LL ww SW Pe Le OWN >» 
<@a ont wes CH ~ jj o— LL LLd wrt Leash © J & je 
pe UL I! FSF EH Wo Wi Uy +e iL mae ODN >u, Il Woz mow ON CONGO se 
LJ @ N ZF =ine D FD ww WWewOs = N<t < ww © LlOttlw-f LL = 
COQ. we _ Amz ZZ O>Qa a ul > WwW rib & et a 
=» OC << Hh mew — iff OZW © wwdnNe OUNN RPwustst Syawo LL 
im & — KC Nw oH — ei te =O be bk 8 ed onl ook AE bed 
— -— et =_ ~ OO» —WZ met LL Om Dey ho m/s LL 
ris OF YW OrFiwO FO WW 2ZOUedO ND RPuS sCOUR LWA Or RL 
>—>—~=090 2 Oper a>O ZO = FZ —BOO AZ NSIS RS Zee SOS rd 
© <—<Iw uJ LJ L) Fp) 
) Ww wn mm eu QO UL. QO b— cima fe eet CO 
aN | A N Ne AN on) aN | NI 
QO OOo ) WO OOO OW oT WOOO Ow 


216 











a JO) OO OCOQDQQOO ODOCOO COO OOOO OO OO O0OO0COO0O0O00o000o000000°0 
PDO PD DS MALI SUN OP DDS MQ PUD OP DW HNO Re Fy PUN OM ONO HAO SF UNO ODO os 
PNiavravia iavfavranroatosTaaraa lan loa tharos log be whe me woot he es oe ew al al allel alialalralialabolonemehenheheheke loin i. 
Va\lallaltalratlalalralalratralral al alfa alatalalral al allalrallal alrallal allel alatralalfal al allatlalral allallallallalrallatial 
d¢ >S>< SK OS OK OS OKO OK OK OK OS OK OK OK OK OK OK OK OK OK KOK OK OK OO OK OK OK OK OE OK OK OK OK’ OC >! OK OK OK OKO OO’ OK 
fe 0.0 0.6.0 48.06.4.04690 646.6.4646.46.6.6.46.9.044.6644.6.6.0.0.46.4646466.64.464.44.4, 
>¢ >< >< OK OK OK OS OKO OK OK OS OK OK OK OK OK OK OK DOK OK Oe! OK OK OK! OK OK’ OK OK OK! OK OC OO OK OK OK OK OK OK OK’ OKO OO OK OK 
DODDODOMNDADMODDMDODMDACADADDDMDDOADADDDDDODODMNDADODODODOOOND 


THROUGH THE 


= 
GO TO 30 


»,BL(I)) 


GO TO 32 


TOWARD 


EOs33. |. 

R Q 

R T 
NT, J 


GO 
E 


- 
GT.O) WRITE (6,52) 


32 IF (LIMT.LTeNLIM) 


Zola, 


A REPLACEMENT FOR VERTEX J. 


RECOMPUTE CENTROID SUMMATION TOG 


BLE TRIAL, 
O) GO 


my « OLLJ 


35 N=1,K 
= SUMC1)*VCLIN) 


DO 
35 SUM(T) 


SUM(T) /FK 








DODDO9DD0D0O0 DODD ODDDD DODD DDDODODDOO ODDO QDOD ODDO DOCOORQO00 
I FAO DO DOSN MFA OE ODO MAM FINOEDDOSNUATINDO HE DRAO AMAIM SIN OM OO 
Re fre eee eM MOO WD ODDONDDARARDAARADHDOO COQ 9 900 O wh tent ett dt tot 
UREA LAA LA AL LAT LA LU EL LAL LALA LALA LALA LAL LAO DO DODDDDDDNDDODDNNUWN 
o< >< >< >< CDK OS OS OK DK OK OK KX OK DK OK OX OK OCD DK OK OK OK DK DK OK OK OK OK OC DK OK OK DK OK OX OK OE KOK OE DEE OK OX 
Q2aaaaaaadaaaqaadaaadaanaaagaaaaananaaaaanaaaananaaaaaa 
oC >< >< >< > >< OS DK OK OK KX OK KO OK > OK > OK OK KK OK OK SO OO! OK DOK OK OK OK DK OO OS OX OOS OK OK OX 
OOOOMONONOOMOD OO MONMOMNNHOOMOMMD OHODHNMOnNnnDMOH MOD HOOOHMMOnooONnNDM 


foes 
Li 
| 
Fae 
tL 
a] 
(a 
LJ 
> —~ 
Pd atAlOm } 
| 
>) Wek uw 
—~— Corr aea 
(Gt se eS | UL LL LL Ld Lk 
aa pant foamed fomng Ot pang 
-~ 
<a rt 
(LJ = bums e 
Ce = < 
o Uj °@ 
=. oe p= 1) e 
ao UU wf 
iL w<t of 
e rs >a OT 
“ < i oo 
ewe WwW VS) p—e 
> cf we ey? 
o O <ax~ — 
~- 2 m4 Lj © < 
Oo => W bn ES 
raz WwW Tp) eYIioO x 
e LU WAC. Lu 
oOo > = O>uWw — 
- 2 WH © Zz WA *e 
~ qc) ~ sWZiUOw 
Oo ww a Lu - OIOOWw> 
Oo OO — co WwW 2oIZzL 
Zz ete = ~ Wee | TOW 
- Om ~ wer = OQ NOray ~ 
zy FO W al eS wad EOL a i 
~ oie also 2 —JtiJ =O ee 
— Ow 2 Ome iu dIiw<D> WN | 
— = © (oad < —- Nt — bh. 
> ee ~~ CeO wW m—™ raatWw lw = 
+ Om~w~ & W © bh re OOO DAW HO 
~ Se 2 aot cr’ ZS Wide ezt5 
ina a ~ Su = $l WiU<r LO <=) 
~~ eZ we DZt > a SS ha esa Lt, 
>= COe 2 2 - Oo ZWitw a tikes a2 
= — —— ~— Oz oe mOOZZuL wwe 
of wa, em 2 YD — CF [seas OP LL woe LL 
= waa =—- we HY ee OS Qk pe > ne) 
oO Wow ww OD AQO-. CO e OO mm © LW => e 
cy ome CA OF S=o+ Ss aMete CVZBVZowanzD>> aoe 
-) -a~ = 40 a ae ra a ae 4s a] niu 
CO dm 2 west OnmO'WA Re etdet= Ww = 
me ee eee eee OLE LOO — oO 
= wa tO il oe Lu z= madd 
=o w. Wwib << Zu i Oa eed CO Qt 
WA} QW 2ySoomem  D w MDD et eee  emmomomom (i) (Lat 
<LL LJ CaO PIN — 
fw © Oo TO 2 ro OJ LL 
Mm cr) vr -P w? ~ 
{o) ©) W WO On Ooo OO OUOWOOON OW 


218 


Ee 
VS) 
26 | 
a 
fan 
> 
Li 
fe 
ae 
LL 
a 
—_ 
z 
<t 
De h 
— +r © 
LJ 
‘ef G2 
Li he WO 
=~ << 
— Oo wW 
Lu co 
fa @} 
om YY) 
>» — a) 
a =z «a 
_ = = 
z -—- cc 
<f +e — 
W a. 
bene BUS Pa 
Ge ~ © 
—_ 4. 
=z uy 
a — — 
_ —_~ = 
Y am fer DCs 
oe mrt LJ 
N~>~> FT = 
—Or ~wo<t WY 
SS Ae 
Za OS a Wh 
OreTtLowyd 
mf) ~~ © =—. 
KvOoz -— D> 
Det Sm 
wah Zz => 
oOo Orexkt~ 
SLL OD Se | —- om 
xc ~~ +O = 
Nem Wt ete = 
wee Oma LL oer 
TM eneZtz TED 
Jem tee LL Le ee er LL 
jerk Ca 
We? 0 OQrweme ZS ii 
wot et QE ew YD TLL 
Te. Se ~ 
pe bee SF Lo I 
ie wee lame fang —— 
Iu, = Pada 
remslLColubin Ws > 
Yo Smee > mw LL 
LL ae — 
“Oo aa) O-sr 
+t +t 
OOW OW 


DO 45 [T=1,NV 





DQDDOOOV OO OOQ OOS OOOO OSS SCV O OO OSC OC 9O9COCOCOCOO C00 0C00008 
DRACO FLD PHD DO HAM PD Ol OPO MQ DPA GS PO DO HOUM FINO PH DDO HAO LN WO fh 
foViav lavratiaraiatiaialafsataslanlad apr araa lad lan leah mh ges homes ome ew ee rd alata alialalralialaliaterekerelenel els 
0000000000 DODDODODQDODDODBDOODODODDOBDDOOODDOOODODOUDOODOODODOODOOOOONO 
>€ >< de OS OM DOE OS DK OK OK OOK OCDE OK OK OK OK OC OK OK et OS OK > OK OK OO OOK OO’ >=’ OE > OE OE OE OK OK OE OE Oo’ >=’ OK >=’ >=’ 
0240.04446.4004,.0440406.446440060060446400040000044046440490040004064 
>< >< DK OK DO OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OK OO OS OK OK OK >= >=’ o> OC OE O>*K > OE OK OK >’ >= > > OO’ 
DMOMODADMAADAODADMNDMNAADAADAADADADNDUNDADMAMDAHDDADNAADDADMODDNDAONADNAOMNA 


NCE sNV eNVTeVeK FUN, V( ly>M)+-1 


oOo 


ae 
fa a 
mm «Lf WY ae 
A fk = WY 
pa octe OS} ae) 
= xr <= 
Cee WO = 
fact = € 
—a oO oo 
“YAW = 
pecs ne ~ 
KOO «< fh 
Pa - Ga 4 LL ~ e 
Lpew f= UY UY 
wor tc = x 
oe ee © eee © LL 
faa ec) Son & 
<T je e © > 
mam tu elL 
“LcioO ak Fa) 
<k= ow <« 
>< fel = 
uy Mec -_ 
Seat LZ 2 ph 
2Zz<ce wet ee e 
™m UW Lh Fr = WN 
~O> KO ~ 
wd LL QOASZnO WwW 
ee Tr 2> ~ 
2A<te “(Oo 
CS J)» We moO 
Oe = AW WY 
Lew pr eOOOe - 
O>— OClzw 
eye 2uwWmN 2 
Zee FS ue O 
O2uWw tik —~ = = 
= <I 00 —Tlu e WO 
O~mam TeKUwNO 24 
wey FYOidIZ2Z >» 
MATL Ww <_J LL. 
mo LU Ho Oe -~ 
GOOow Le IM eeniv = 
‘| my eT fh S> ~ 
=> Cc) x< <em 2 
ZeZ YsmnZzZzwnNnke D 
<[rtewrt Tet Ohm et) 
OU, st wet ey Wl ~ 
>< o> A od eed BO on A So 
Uitte Fase O ~O 
C6. 6) eos © tA 
mt FCN Sr OS 
Oooat ocscouw m— Ws 
905 2 SCS) a af ae ae af oo LL. 
C00 OFOrO0OrFO UI 
Wine «= F-Ne Ze 2a 
want Meal Sag? we ee set m— (_) 
LL — =— 
ere per lee pee ke bet “WY 
Qeiet Arad sest OZMe 
= =OFtsSz= = Sw 


—_—=—= 
CAeeaCIeKer COO ii 
90020 OQ00e 022— 


Mee *LOWILuN LWMnost 


—> nd 


COm N MPN O 
FIM MW WWW WW 


PALES) 


GO TO l 


).LE-FUN(I)) 


VeNVT eV eke FNeCe IK) 


=1,NV) 
NVPyNVT) 


(JeTded 
1 


J 
G 
(V0J,TdeJ 


GO TO 2 


TF (IK.NE.O) 





a —— 


DOQDOODDDDCO90N090009009090 
DDO ASANO PAOM DDO AUN FIO 
DOR Bette Rothe tee heh OW DOW OO 
DOD 000000000000 0000 
SSL KS DKK KK DK DK OK OK DK DK DKK OK OK 
nooadaaaanaagaaaaaad 
SES DS HK CDK DE HK DK HK DK DK DK DK DK PKK OK IK 
DoooooOoHnnnHonoonoonaea 
a} 


Q 

—a 

(lu 

2 
Ped LLS- 
e <O 

> 
UI) 
Sm oad 
wz oO 
<{-—-+<f 
pant ef tt ily 
(a Mo a 4 -~ 
pe poe <E ~ 
n> ® 
Ls Zz + 
OF — 
nal SP aw LU 
tang LL) am fo 
MM «OQ ~~ & 
aqOz am eX 
Uj af = 
Le free ond (\} 
mL o LL) ~ 
xo = TRk™ 
orZz e a 
Zee M UIKXK « 
e ie r [~~ ont f° 
eum @ ~f PAl Hd 
x eS mt wl 
Pe OQ Ww ash 
ee eo Cif @& 
Ai= ~~ we OX 
— Ww x SAN 
n> = ea & 
eZi NN elu 


O<{emmwws Sh 
ji=—>Se NN Ste & 
hm Wh ilxX<e 
Met Zhe @ hh xX< 
ae pe | on Pe» «> 
me pe <fitj mix) 
C>OoOtr swWwM 
mujhe o& he LU 
2”Nx Act 

awd Zl Nett 
aD OO eO>r0a 
poet CM Cl mm 


=1,NV) 


TO 3 
I) ,1=l,NV) 


(629) IK,(CC1) sT=la NVI 


O mNaTteOOO 
e ON Set AQeia © 


—e < Or preed wast a Tae et Sage? 


(697) (CCT) 1 
( 


- = 


fea 
C4 WILL OF LU Od UW Sinz a_oooo0os 
xe ex 5° Se ta Le ees LU he LL Lh, bd 
at AICO 
N an! fP LN Orc OF 


220 





FUNCTION DEGRAD 


This function was programmed to ccnvert degrees to 
radians and radians to degrees. A third purpose is to 
ccnvert degrees to a range of (0-360. [It is used extensively 
thrcughout tke programs, functions, and subroutines’ listed 


in this thesis. 


221 





FUNCTION DEGRAD 


Ne (@i@ 2 @ leo @leolelel@els)lelele!@(elelelelel@ele(®lele|®\eolelsiea) 
MAI FIN OP ODO SENHA SIN OP OAS AIF PFOOM ADOC 
oad goed ad peed ed ond ered et teed CNY NI CUCU CU CU OUI QUOI CI OO 6 OD 


AD AQAQDADODODOHAOQOAOAAOOAOSOOGOO0000000000 
(e Mle Mle Mle Mla Me Ae La dle Ao de An A> Se do a Me Aa in dle Ae Ae Ao a Mia Ae a Ae Ao a Ao ao 4 
OHNO OOO OOOO OO OOO ONO OYVQOOOVOVOOUNOOUUNO 
BAAADODAGDGAGBOASGOAOOAOGOAAA ARAB OAAGBAAAAAGa 


© 
e& @ 
Oe © 
—& © 
S&S © 
WO 
=a a 
{Zz UJ 
m= LLC LL) 
Qujs<t = 
<{ > —J eal 
OC he. LW 
Ww co 
Amd 
eae UJ 
<I UL) = ©) 
co = =. 
° MZO <f 
m 2 fa 
= <r J 
QO — OO alee) Lu 
za O<ctiu ee aa Be 
=) <r LL <f — 
bw mwOuw _ 
~ je OO ‘on +r Oo 
Oo Oo = W<f Ue 
—a  eY yee Oo & 
LL tu = At @ 
rc NWeE Oo a) OO =m a 
SY Uc —— Lu mis 
az Wor (am) — -F & uw 
- aww WY _ Nett WW Me 
oo OW awk WY — WO © 
=e tk 6t) qc WwW ¢ S@Qre Tn! Ge — 
a Ow cmv tus heey ON OD 
ce re Qa0o> 3oO> ~ WN on) 
CO Kee Twsrz a2 eZ O00 =» OF0 OO 
=_ ec ~OoOO <ooQ ede SF OO e ® 
— Whe (2) OW Oo = Oo —O 
> (LW om fame fmm OF aw ome CO) eo mO0O0 OVO 
CO weawy Lot Ll lk Cie MO at OM Xan) 
=< Ord wa woo eet +¢ mi or D 
eC wom >>p>2Zz S>azz FOOD DPD GsD c=) () 
5 ae za OW OO Weo WoO 2 
WwW OCOO=mM Ooo O02 O Z2Zee2tz 2 ee Pl a =) 
Swe = FOO OW GO DOOD DPD YO eD> o> WwW 
Ioo eat Ss Su. | 
me ome = =6f | j | LL. LL. LL.) cS i 
Oo Qnore WECUNVHEODSZUIMNZW S 
rt LI Oe OO amet COC met COQ A> ae ats Pa 
— RUM es SaeaZaODaO2il [Dow lode 
OO Oct mit LL CL) weer ase) be CD pre see see) pee eer OO = 
= Zooawe ake = ee ae Pal =a OfKO 
=a It we WY PIULLOODWDOLWLOIOWDOWW Ss 
tet wodctoaz Zz Ln pment oo LL LLL CO et toe LL ee LL COC LL 


sae CO —fP uN 


WOOQYOOQVOVUOVDOVUO 


Bee 





FUNCTION DELAY 


This function is used as the time delay in the speed 
ccntrcl cptimization runs. It was designed to be used as an 
equivalence to DELY in DSL Simulation.The following 
variables are defined: 


E (1) is the storage array (should be initializeéd 
before the first functicn call) 


K is the delay step count 
SPDLCES is the variable to te delayed for K sters 
P is the flag for delay or no delay 


P>0.0 delay SPDDES 
P<O.0 function cutput equal to SPDDES 


The functicn stores the input value (SPDDES) in E(M) and 
decrements tke value in array E(I) at each call of the 
function until the value is in the positicn of £(1). The 


value is then cutput from the function delayed K intervals. 
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FUNCTIQN DELAY 
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FUNCTION FE ~- RUN A 
FEA 


This function is the simulation for heading ccntrcl 
Optimizaticn of the approach phase. It is called by 
subroutine EFOXFLX. The integration step size is 0.04 with a 
final tine cf 20.0. In this function all initial conditions 
are set to zero except initial geographic locaticn and 
speed. The reference ship maintains a straight course and 
the contrcl ship starts its approach 5 shif lengths astern 
and 0.4 ship lengths laterally displaced to starboard cf the 


reference shir. 


The function is referred to as functicn FEA i1n the 
text. | 


Duh. 
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FUNCTION FE - RUN B 
FEB 


This function is the simulation for heading ccecntrol 
Optimizaticn of the turn phase. It is called by subrcutine 
BOXPLX. The integration step size is 0.04 with a final time 
of 20.0. In this function, the following initial conditions 
are non-Zerc;: 


contrcl ship rudder angle D2D & ¥(14) = 8.7 degrees 

lateral displacement Y(11) = 0.2 

reference Ship's speed 001 & Y(4) = 1.0 

contrcl ship's speed 002 & CDCT2 = 1.5 (after first 
step beccmes 1.0) 


The reference ship's rudder is activated tc 5.0 degrees 
between time 4.0 and 5.0. The runs were for port side 
replenishsrent. 


The function is referred to as function FEB in the 
text. 
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FUNCTICN FE - RUN C 
zoe 


This function is the simplified simuliation fcr speed 
control crtimizaticn of the switching functicn SW. It is 
called Ey subroutine BOXPLX. The function shown is for 
approach speed of 1.1 and a replenishment speed of 1.0. The 
runs were made for various realistic combinations to cbtain 


an cptimur stitching curve. 
The rcur used a step size of 0.04 and a final time of 
10.0. The twce ships were run linearily with only tne 


lengitudinal direction and motion of any ccncern. 


This furction is referred to as functicn FEC in the 
text. 
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SUBROUTINE HDGRAS 


This sukrcutine was programmed to calculate the desired 
heading (FSIDES) for RAS heading control. It uses the 
outputs of sukroutine RBMEAS to calculate this heading with 
gains RSENS and WISENS. The large number of outputs in the 
subroutine call statement were made for ease of DSL printed 
output fcr tracking of simulation accuracy. 


The sukroutine also incorporates a loop to avoid 
computer precision problems in the ARSIN function. 
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SUBROUTINE HDGRAS 
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FONCTION KE 


This function is required by all optimization runs. [t 
is the function that contains constraints for subrcutine 
ROXPLX. No ccnstraints are present, conseguently function 


KE is set egual to 0. 
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MAIN PROGRAM FOR FUNCTION MINIMIZATIONS 
MINIBX PX 


This is a generalized program which calls sukrcutine 
BOXELX. It¢ main purpose is infut and output of the values 
reguired in the optimization runs. This is referred to as 
MINIEXPX in the text. 
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MAIN PROGRAM FOR FUNCTION MINIMIZATIONS 
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SUBROUTINE RBMEAS 


This sukreutine measures the range and bearing of the 
ferward and after stations which is reguired of sukrcutine 
BCOGRAS. This is done with trigncmetric functions as shown 
in chapter II. The subroutine is specifically designed to 
circumvent any ambiguities usually associated with these 


furcticns. 


It is the basis of the decoupling of the two RAS ships 
that this thesis is based. 
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FUNCTION RKLDEQ 


This function is the Runge~Kutta-Gill forth-order 
integraticn used in ail optimization runs. It is prcegrammed 
locally and is part of the IBM 360 SSP liktrary. A full 
explaination and description is shown in the first few pages 


of the functicn listing. 
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SUBROUTINE SLOPES 


This subroutine contains the table look-up and 
interpolaticrn scheme for the interactive forces and mcments 
presented in the RAS environment. It is long and must be 
pre~compiled fer mest of the DSL Simulation programs shown 


in this thesis. 
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PUNeGTLON (Sern z 


This furction was designed to aid in initialization 
problems asscciated with the DSL function DELY. The effect 
is that the functicn initializes the delay lcop until it can 
ke self-suppertive. 
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FUNCTION SPINIT 
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FUNCTION SPDCTR 


This function is the speed control used during heading 
control develcpment. It is used directly as the speed of 
the contrcel ship with information cf the speed of the two 
ships and the Ilcngitudinal position ADX. It contains a 
linear furcticn at +1.0 ship lengths to a dead zone of 
40.001 centered about the alongside position (0.0). 
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FUNCTION SPDCTR 
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FUNCTION SPDOFC 


This function is identical to SPDREC except that the 
ability tc offset the alcngside position (0.C) is 
inccrporated. This is the speed ccntrol function in its 
final develcfaent fcorn. 


Z15 





FUNCTION SPDOFC 
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FUNCTION SPDREC 


This function 4s similar to SPDCTR except that a 
switching function is incorporated. This is the furction 
used for cptizization of the switching function and 1s used 
in the velccity loop simulated in the velocity ccntrol 
section cf ckapter III. 
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FUNCTION SPDREC 
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FUNCTION SWCL 


This functicn contains the fifth crder felynomial curve 
fit fcr the cptimal switching positicn of the speed ccntrol 
COP. Its range of values for SPDDIF are 0.1 to 1.0 
nermalized speed difference between the two ships. 
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FUNCTION SWCL 
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SUBROUTINE SWITCH 


This stkroutine contains the gains and méchanisms 
reguired for the adaptive gain schedule developed in this 
thesis. It includes the optimal gains cbhtained frem the 
heading ccntrcl Optimization runs. 
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SUBROUTINE SWTCH 
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SUBROUTINE SWTCHF 


This sukroutine is identical to SWICH except that the 
turn phase gain VFBG is relaxed to allow for cffset 
longitudinal position placement. This is the adaptive gain 
schedule in its final form. 
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SUBROUTINE SWTCHF 
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SUBROUTINE TRANS 


This stbroutine takes the lateral and longitudinal 
geographic displacements and converts’ then to actual 
displacements referenced to the contrcel shifts head. This 
is done tc gain a more realistic reference for sukrcutine 
KREMEAS, sukrcutine SLOPES, function SPDCTR, function SEDCFC, 
and function SEDREC. 
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SUBROUTINE TRANS 
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FUNCTION XLIMIT 


This function was developed to allow the LIMIT function 
of DSL tc be incorforated in the optimization runs. [It is a 
Saturation amplifier with a gain of 1.0, and upper limit of 
UL, and a lower limit of XLL. 
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FUNCTION XLIMIT 
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APPENDIX B 


The final fcrm of the simulation program, with all its 
Subroutines and functions, is a very complex and complicated 
maze. Tc aid in following its progression, this aprfendix 
ccrtains a detailed block diagram of the program with eéach 
variable listed in its computer variable name. Each page 
contains a functional part of the simulation with inputs and 
outputs shcwn cross reterenced to their origin and 


destinaticn. 


The fcellowing is a list cf the blceck diagrams ccntained 
in this appendix in the order in which they appear: 


Ship A (Beference Ship) Simulation 
Ship A Heading Simulation 
Ship A Speed Simulation 
Ship B {Ccntrol Ship) Simulation 
Ship B Heading Simulation 
Ship B Speed Simulaticn 
Subrcutine RBMEAS 
Range Measurement 
Eé€éaring Measurement 
Subrcutine HDGRAS 
Heading Cortrol Loop 
Auxiliary Functions 
Yaw Conversion 
Coordinate Conversion 
Feedback Loop 
Rudder Y*cdeling 
Ship A Rudder 
Ship B Rudder 
Wave Generator 
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Wave Direction 
Rardem Variable 
Wave Enccunter 
Wave Components 
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APPENDIX C 


Thrcughcut this thesis subroutine SLOPES has been used 
to output the interactive forces and moments between ships 
in the RAS situaticn. This subroutine, adapfted from ref. 
11, does net cecntain a complete picture of the circumstances 
envisioned. In particular, ship's speeds cther than the 15 
kt. operating point and different ship lengths are not 


accounted for. 


As stated in-chapter II, the speed modificaticn factor 
can easily Le applied for both ships at the same speed and 
otker than 15 kts. with the following expression: 


2 
Seur = CDOT 


Ships replerishing with different lengths can also be 
inccrporated aS shcwn in ref. 1. 


Subrcutine FAMIC listed in this appendix inccrrcrates 
these twce ideas along with a better method cf determining 
the interactive forces and moments. The curves of figures 
II-11 and II-1Z2 were quantized every 50 feet of DX for ail 
the DY curves shown. These pcints were then used in the 
NPGS XDS-¢30C digital computer and AGT-10 graphics terminal 
tc obtain a family of best fit curves. The kest fit 
criteria is tased cn the sum cf the error squared at each 
gGuantized feint {mcdified somewhat by this researcher's 
€valuatior cf test fit between pcints to eliminate spikes 
and cther ancmzalies). The results of this curve fit process 
is summarized in tables cC-1 and C-2, which includes 
takulaticn cf the best fit criteria. These pclynomial 


coefficients are based on the DX distance and are ccdéed in 
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subroutine FAMIC as YY(1) thru YY(11) and YN(1) thru YN(11). 
An interpclaticn algorithm is used to determine the forces 
and moments at DY points between the curves of each family. 
Although all the computations are based on measurements from 
the ccntrel ship (ship #2), the interactive forces and 


moments are also computed for the reference ship (ship #1). 


Figures C~-1 and cC-2 are the interactive forces and 
mcments cutput to show comparison to figures fI1I-11 and 
II-12. The speed of this run was the operating point of 15 
kts. The ships are of equal length (527.8 feet). 


Linear interpolation of the interactive curves for 
greater than 150 feet DY distance is accomplished from this 
150 foot curve te a value of 0.0 at 200 feet, It therefore 
assumes nce fcrce and moment are present outside the 200 foot 
range. All forces and moments for DY distance of less than 
50 feet are taken as that of the 50 foct curve. These two 
endpcints are by nc means exact, but will suffice until more 
detailed cata can be gathered. Another inexact endfoint is 
prceduced at the curve families limits of +550 feet. At 
these points, the forces and mcments are forced to 0.0 since 
detailed data outside of these limits was not available. A 
Side effect cf this abrupt truncation will manifest itself 
in the instantanecus commencement of the fcrces and mcments 
during the arpproach phase run. The endpoint variations in 
scme of the curves of figures C-1 and C-2 are due in part to 
the curve fitting routine used, but mostly to the 
differences in computer precision. (curve fits were 
calculated cr a 11 digit precision XDS 9300 while the curves 
were plotted cn single precision 7 digit IBM 360/67) 


As previcusly mentioned, the speed mcdificaticr for 
Other than the operating point cf 15 kts. is cnly completely 
valid fcr tke situation where both ships are at the same 


speeds. Sirce this thesis considered an approach phase 
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Curve Fitted Interactive N Moments 
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where the control ship enters the interactive field at af 
speed guite different than the reference ship, some 
modificaticn cf .- the interactive effects shculd be 
considered. However, exact relationships are not available 


to compute tke required modificaticn factors. 


To dispel any problems with the design of the heading 
ccntrol system, the worst case speed modification factcr was 
chosen. This factor, in effect, considers that the 
interactive fcerces and moments are derived from the ccntrol 
ship. This is accomplished in subrcutine FAMIC with the 


fcllcwing fcrtran expression: 


SPLE2 = CDOT2**2 


As stated in chapter II, it is feit that it is more 
accurate tc consider the interactive forces and moments to 
ke modified Ly the speed of the reference ship, and can be 
coded in svukrcutine FAMIC as: 


SFLIF2Z = CDOT1**2 


With the scenaric followed throughout this thesis, this 
expressicr wculd equate to unity throughcut the RAS 
Situation, since the reference ship is maintained at 1.0 


ncrmalized speed (15 kts.). 


For the sake cf error analysis, Simulation of the werst 
case modificéticn is performed. This gives rise tc forces 
and moments 2.25 times what they were in the rest of this 
thesis during a portion of the approach fhase when the 
ncrmalized speed of the control Ship is 1.5. Figures C-3 
and C-4 show the interactive forces and mcments for the 
approach rhase cf the simulation. The comparison plots 
which appear in chapter III as figures III-24 and I1i1I-25 
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illustrate the extent of the changes. Most notable is the 
smccther CUtERS of subroutine FAMIC. This more 
realistically portrays the interactive effects in the RAS 
envircnment. Figure cC-5 portrays the geécgraphical flot 
which ccmfpares with figure IITI-26 without speed 
modification. Although differences exist, figure C-5 
illustrates that the interactive effects speed modification 
factor fcr tke worst case does not drastically alter the 
approach phase cutcome. The heading contrcel system cesign 
is still valid in the face of these changes. 


For reference, figures C-6 and C-7 show the interactive 
forces and gcments in the turn phase as calculated by 
sukroutine fFAMIC. Figure cC-8 is the turn phase lateral 
distance flict preduced. It can be seen from this 
illustration that the maximum excursion is 0.0056 normalized 
distance (2.S6 feet), well within acceptable limits. 


| In sunpmary, the designed control system will accomodate 
even the worst case modificaticn of the interactive effects. 
This insercsitivity to a large range of perturbations, makes 
this contrel systeg a more viable design fer actual ship 
installaticn. 
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COMPUTER PROGRAM #1 


This pregram incorporates the ship dynamics of two 
identical Mariner hulis. These hulls are superimposed in 
space to ellcw fcr comparison of the effects contributed to 
rudder mcdeling differences. [In this particular run a step 


and ramp rudder were compared in chapter II. 


Ancther kenifit of this pregram is tc set up the two 
identical ships required for the RAS simulations in chapter 
III. Basicaliy, verification cf the models in three degrees 


cf freede€az is acccmplished for the Mariner hull choSen. 


The glcts produced in this run are shown in figures 
mi-2 and IIi-:. 
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COMPUTER PROGRAM #2 


This pxrcgram models a practical rudder response fcr a 
Mariner shir type. The rudder limits (stops) are set at +30 
degrees anc the rate of response is limited to +2 
degrees/seccnd. A scale factor (LUC) is introduced to 
modify the respcnse to match real time of the mariner hull 


chosen. 


Twelve passes thru the program are accomplished to 
cconfcrm to different sets of initial conditions and final 
desired rudcer conditions. The plots produced in this run 


are shown in figures II-5 and II-6. 
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COMPUTER PROGRAM #3 


This pregram models a reduced order (first order) gas 
turbine froprilsion plant for an input-output relatianship. 
The program does not scale the plant to the mariner hull 
used. This was done when introduced into the main 


Simulation prcegram first listed as computer frogram #8. 


The time delay (P) is assisted in initialization by a 
dual feed into the system; one thru the delay itself and one 
directly intc SPDIN. The programs can be modified to ccmpare 
a family cf curves. by introducing the following sequence 
into the TERMINAL region: 

INTIGCER NUMB 

IF (NCUR.EQ.NUMB) CALL ENDEW(NPLOT) 

IF (NCUR.NE.NUME) CALL RERUN 

NCUR = NCUR + 1 
where NUME is the number of curves desired (less than or 
egual tc 10) which is set with a PARAM statement. The 
comparison is dore on the conditions set in the terminal 
region [i.e. decrement or increment the system gains (€g. G 
= G + Q.CZ) }. 


The plet procuced by this run is shcwn as part of 
figure II~-10. 
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CCMHPUTER FROGRAM #4 


This pregram models a simplified wave Simulation 
ccapesed cf two superimposed sinuscids (fundamental and 
seccnd harmcnic) and a small random wave. The aodel is 
inherently scaled -. to the Mariner nondimensional 
characteristics. Introduction of these waves is accomplished 
in computer program #7. Subroutine DEGRAD is shcwn in 


appendix A. 
The sea state force plots in the dimensions cf the 


three degrees cf freedcm produced by this run is shown in 
figures [I-17 thru II-22. 
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COMPUTER PROGRAN #5 


This ficgram uses the mariner hull model Paes c 
introduced in computer program #1 and the control system 
designed in chapter fII to simulate the afrroach phase of 
RAS. The subroutines and functions that are to be inserted 
frcem appendix A can also be dore in object ccde by changing 
the word FORTRAN to OBJECT and placing pre-compiled decks in 
the same locations. In fact, due to the lcng length of 
subroutine SLOPES, this must be done to be able to run the 
Sipulation with the DSL default job control language (JCL) 
presently irstalled at the Naval Postgraduate Schocl IBM 
BioiG /67. 


The flicts produced by this run are shown in figures 
IiI-7 thru I1I-13. By changing the gains and introducing 
the foilcwirg code, the plots of figures III-14 thru III-19 
are prodtced: 

DILES = S.0*STEP (8.0)-5.0*STEP (9.0) 
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COMPUTER PROGRAM #6 


This pregram combines the approach and turn phases of 
ccmputer pregram #5. The added subroutine is a result of 


Simulaticn requirements to switch between adaptive gains. 


This rrr froduced the plots of figures III-22 thru 
IlI-34. By substituting the initial conditions of table 
IIi-3, tkis program produced the plcts cf figures III-35 
oad LItT-¢4, 


Sy ere! 
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COMPUTER FROGRAM #7 


This pregram combines the calm sea simulaticn of 
ccmputer fregram #€ with the wave simulaticn of computer 
pregram #4 tc Simulate the model and control system in a sea 
state. The waves are introduced thru the rudder 
nendimensicnalized coefficients as shown in chapter II. 


fhe flets produced are shown in figures ITI-66 thru 


III-73. Figure {11-65 was produced with the same program by 
setting wWl=1.5. 
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COMPUTER PROGRAM #8 


This pregram incorporates a fifth crder pcelynomial 
curve fit speed ccntrol switching functicn te give optimal 
lengitudinal fositioning. The scenario is the same that was 
used in tke cesign of the heading control development. The 


low order secdel of the gas turbine propulsion flant was 
used. 


The flets fproduced are shown in figures III-8C thru 
ti1-383. 
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COMPUTER PROGRAM #9 


This pregram introduces a longitudinal pesiticn cffset 
capability. The sethod takes the control ship tc the 
alengside gesition until 450 seconds into the run. After 
that time, with the ship steadied, the cffset position 
desired (XCFSD) 1s switched to the desired offset. This 
method negates some of the transient oscillations which 
cause urstakle conditions in the apfrcach phase. A 
seccndary change is the use of subroutine SWTECHF instead of 
Sukroutine ShkTCH developed in the heading control section. 
This new subrceutine relaxes the heading velocity feedback 
gain (VFEG) tc allcw turn stability in the turn phase. 


The flots produced by this pregram are shown in figures 
mer—e8 thru Iiz-101. 
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COMPUTER FROGRAM #10 


This pregram incorporates the sea state ELESe 
prcgrammed in computer programs #4 and #7. The WX wave 
force, hcwever, is introduced at the end of the propulsion 
locp: to allew more realistic perturbations. This is the 
final fcrm cf the complete heading and speed _ ccentrol 


systems. Tc cun this without a sea state, set WFMA tc 0.0. 


This pregram produced the plots shcewn in figures 
mets thru IIl-105. 


562 


COMPUTER PROGRAM #10 


Iz oct 
Oagadoe wn 
DP) od ZINE < 
Te) ——acCq | or 
T — 
LL LY Ow. ha 
= 1° 8) eOUORLO 
—_  ) QO awe 
i © OO wm 
a ny N2AaZ > 
- oO LaAO~m 
wy it Zz Am< 
ff an) Or <ake=zO 
ve) pa —O miu 
— x @& Fe l(LiLY) TS 
~s ® — OOO 
© Le Oe & NY Aegr_) 
= OO ACO Il OOwW wy 
WV) —j OOonc oa. JO} 
bh Ooo Wu dou 
= TS) oOoodw < NO 
_ I< eCOO0Q T Lev e 
cw fae 4 @eex ® OWWnNOO 
= i LL UL 1oOo e (a) WO > lw N 
=) _—_ wu two iT N= ODS Ov 
= Zaz COUN NI 7 NY OW a | : 
o — pne >xvae © Oo JOoOOW<<+ e 
-~ NZZLSO > bh @@QZLuLzToo © 
oi Lis LL em Nl Ke - N 2O Ot i 
+r >> ZmHtms = @ ro) Z CIQMN20 C5 
<{ <<{ WL Ne ¢ ~~ Odw alwW 
LLy xk 4 Om ae | © - Zz Noe 
o~ COOWNO II " 2 0 ae ee 
Ww 14 Oa OQOoocro sn SoS Q0 miWw 
~ MUL eO ey eZOM tt W Wena 
~ =poJ WLS QwWreOxZ2 So YCZMOF DO. 
© OO NwWwIi OGgiaqaie«+O 2 sag ited ~ 
~ ww AOw wee ReOm KNWWBZNOW © 
~ —— N mO>aD sweet ek © elpi<c re®) a 
oO ae NO YW 2K Corfe, eHONO > O en 
~ Oo POMNERO eS envio Gh ae doz NI Ea 
N Cc) Nom Pwo cel eMOOAntZ ww ww Det it ee) 
— «6 OO cOMDHOZNNOLTS WOO NOrtrteuzZzoo— <r Se) 
QO 0D esti OINAIN (OM>~CO>-D>rFEOOOZ- OC OF Oo WW 
M Owl ail eNPONZOOOS e er LK Leo em e OO 
O Albw N tO! «€ eM>~O eO!1NIU]OIUO COM mY (Or WwW ee 
= A0 & lH COOUMNONST HO eH ell I< o<z ii © emt) | tee ZO ll eOnQoee- 


ANN KORY CNZNIHOONHOCZ SOK HON WwW &ECLTOVOs oO eli t 
OnMD LdOm a I OO | QO!) YOR NEO INZWY Ee ZOUNNNN | << Il 
LOAM ws HWY ZOrO Il > hee Se Sete of 
me 2a VY ZZ NNR OLN Ye ~ POU ZO2Zza>~OoO0 O 
ZOKweaee BZ SZeCSCSoR IZ ZT > ORM Xe >> DK 2 


m~w!lL) @ So ae 

COQh WOW RE WRESS eee eee et 
Cw YusOnowndaatatdoqt Nnnnw 
D2 OO BwmeRRZE2eeeKerere 22227272 
~SN Sz CO ZO ttt ONTMO0O0 


a 2.22 
gd oon 
OQ YOVO 
2 222 


LiL HW a ON 
OOHNOUIKeWY 
<xKISFSSSS 


SSSTsTTS {T= 
Saodadaqaqa 
Indo Ao Me Me AeA Mod 
aqidddqaq<(<q 


SOS Fee tO et 0.0200 OO HOO OO F Mt mete He 36 Ht He He HHO OOOO 


363 


0 
X 
*A2Z1 
5 
1 
1,1,02D) 


DO1/K 


Oonwd (aa tromnnrea ~~ >) 

> < ZO. @ AYCYN ena OO VDowna 
{'{<«<OfF OA A ke AO WN er O= ex O Ue e/-~ 
NNY I DDS eer > rZeDDOWWAPpoRsn *eQOO@w~e Nie jj 
OO RN) IS Oh Oe ee ZN ZC KX HOOF NO COM 1OOAOH Ie 
~~ DPiriBr ir St ta mnt SVs t SxtKaTVaqnysow ti ~woeaNiy 


LL 

= 

3 

© 

6 

WN 

=< 

aa] —_ 

= NM 

~ 3 

Ca) et 

=, ey, ah 
© Qo = 
® => ™~ 
© Lis N 
4 = + 
4. WY) + 
= - as) > 
<I © W) LL. 
in ® + = 
+ © > + 
—~ ox a Se) 
-~- ™~ + N 
Ont ian ON 
om = wy 
OuwWL © = 
Tz ‘am wt 
“<I O On ~~ ® 
m= e@ O+ (om eo 8) 
aLS =O = = 
+s o& WIA ii + 
e ~f\ eQ ee SE ee oe 
Och AReN OOC0OW 
we OO wet SF ess 
Nibhet Ost WA 


=1u4snodta eszZzZzy— 
LL wer <E w <I (1) bs pee ss CY) 
See I CHE NNOS 
wert OL A OOO OD St Ht HE HE 

=O. Sus etLIL UL UL iL 
to 2) NH ONWSSSzZ 


Ct IU HEHE Nn ban ae EORONL eOCOw wht FON Sat 
OO 434 ANS Stat SCN NN NOOO Seat << L) Hh OH ODTOOO Spe Hi > HOOMK I Yr Z KILL 
$A TT QOS ANNES ANNO DOO HWW OMDW ANS FP NOAW et LLY > SSW LiL > 

AOQOOVUVVSZaodtaadtoaat odo yYYOOQCOUYOOOYANWYHwR SSSFWWSESSTs=z 

= $ 

— _— 

— fad 

ra Lis 

_ +f S 


364 


LLL 

z= 

+ 3% 

oo 

@ 8 

MNIAI ad 
a eee : 
Pile a ad 
id om ‘:l 
MY) ww © 
% 4% LL © 
ae eee | NX <x WY 
zz a. —_ om | 
~~ — =—— Te A 
NIN VY? ~> wmac~ 
se = ~ mM co alo) 
 % — — an oe NOU 
zx< om UJ) zZzMN QOOQ 
UL UL ii f= -_—— me () NN COaO 
ss HY wot aN WW) mt LAN 
te +e Mus O+ LLL % + LE Le w!wn% 
00 + CE oa meme COC) | — 4 =) CO OL we 
NN Om Siw -~+ “SN bk b— ++-~N =~ 
Or oOr eS tua oe -_ OO NN Onn 
LVL We ee ee tN mame OO Reem OWS 
sui | © J “O00 Net ast Se ~ZNOONANN CO 
PP PiUlwtue «© MO OOO at eeEO | + ZZOOOAOKMTeEREeES NOK 
tee (SS ET IOS ott 2zauo ** OOOo-;-;--— HHZOOXKXOC—™ A & 
© ewww A OT 4) +HtHeteats OOO aaa mete HH OUNOONKWY 
NONZ2ZZe © WY aaNONet OOUOKR rR eE LO AnnNsOooonw 
—u NIC) fF «2 < SCUN ile <[O eQOOodooo Orv YY NA AIaSe el OO < «OD 
+r RH «€ JO Onm= |! | e eet Q—~— OO -++#QOM@Md! I « OA~0aC 
ome tt HE OOO =O ee — 1! | ees s *O0on2zx<> ALANA | | Newt ee OO Wf 
WI ZK OMe! OC oe — QNasaMtUIOOOD Oe & = ZINNNNSROS ALhyh—~s 
SSiILUiLw | OW et Be CO fw fem fe ef <I ewer DY) tet alo OY Od Oe 
www SSS 1 EO Fp) GTOOOF# #& PUI UOt + OO N- $I Oot %H UICO>- 
ZAt He HN Re Stew et steat tt OOF tet OY OO wr eet OS WINH QONN GY COA <{e4 
mS > HSS ZOO aR BN SOOO SRE eww ZOO egqaNoOooans2a 


NA LL CLOZRea wo OC) Ht H% Ot me MO Ee EL OO SI mO te 3 Ue it et YL es 
Se SS et J OD ed tes ett YY SIN) wrt FO ZOO eet et YN were FO 1 O00 
ZK eH PINON RE fH KI EE eee OO WR KO HW Le 
LLibibibin gp Wt 1 I ISMN OOO See HU RR Oh eS I OONN YY EONS 11 
Se KH SS ee I 0 ee eS eset et et 2 Nt to Pi emReNNSww 
be SSSWNWO | OxXKxye Dasa i i OCOR RR le Ree at > ODNNN IH HOOrRFMNOOWw 
Wih ww hOR RRO FH ES ANGN SSVI ABOQOOWUOOBOO ! ' SAICSVAnMNVNNOOOO ! Ong 
ZK 7K wi OC tt ere IL IL Let MOO QQOWO KPO ORR SDRMLLIiLaSNOOdOuw £0008 
SS FSFDBOOO OOOO New et ot TO <f OL) OK KE NNN et es et ee <I OO KT OWI) 1) 


365 





ERR*¥LUC) +KS2*WX 


LU =) pe om me 
ON eONNAIA 
ORFNO> +e e 
AOOnaaadoO 
OOD ~~OO 
+-+$OrrNNN ZX > 
ai JOC * « 
WZ KPOMINA 
OmODH HOO 
Oo OeRKNALX< > 
MA bk wwe 
—Ne COO 


N 


%e— OOQOCAANX 
SH WOOO O > > Il 
ne 2itke—ROoHn 
ZANMHNN SZ WNW 
tie 1 bem ee A = 
Oocnoo i #sae 
AAO KQOONNA SD 
MM OK Se Oe NY) 


NOSORT 


am pee ome (NJ ord CO 
aN eRe MO 
Sel SQND&M -OO 
I~<iteod esdsom + —-} 
we ce O ZA AKOD eDO-OoOMmM 
wot oe Set Se ee @¢ & 
OOH et CNO ee O OO 
HS -O~W>> a HOOK OW eo 
Ow Ow 1 © Iwo ONDHO 


OO e« 


~Ownwn ONZS2Z ODO 1 M01 Ne 


O<OdaC~~ack | 


ee! N~D ea 


Ivy avCoOaNnNnoconkr~OkF-~ = 
COCOBNT<ST CONDOO UW |e ~ 
OmOw fw RW Se ens eo 
WOWO-YZzOSO>IOODZ eR Ost 


OCOvwO nw i 


DO tt treSTrtetkee 


Hen NOweTlaaga 7 dDD yO ZO 


soo ol 


WN OOO it tt Wel 


OOKVOO004I0Or RK SOOMWA 11 O 
SSSSSOYNUMOO DW we EEO tt 
Ieee INOOOWM YS IN 
>> COCOA OOO SS OOA0O 


366 


SORT 


»DISTES/25.0) 


EPARATION 


Mor ao 


>< Lo 
M—Oni 
A N<t 
mer 
TY <1 
YM «or 
rt <I fr 
> LUI 0) 
Ow QAn— 
I SNZVuUWO 
mp De UI 
<W>S «© WY 
OM -v~ 
<< aro 
~NQefezw 
~zrh-O 
Ome ort <le 
HH - OLW) 
KO e<etm 
Oqt<toOaw © 
BS ~ etli> @ 
seHNx<eunz A 
SwO 3+ Own 
IOdr et Ola=™ 
Smee te tN 
we (UY) om Ht LL 
AULWOF SH M 
=. O <f m—c)+ 


Eizo KKK TtiMn~]~ © 


WN | 


i= 2NnOore—S I oO 


et HOOKING G( At l~ 
=H KNW < ew Se SSE il 
Mm He WN oe Oe ee 
—s—OeM OX at heILOOrA 
Oogooodt Oo OwOmsLaI << 


DYNAMI 


OQ 
om! 


+ 











eWTSENS,RGN,VFBG,BDOT20;, XOFS) 


A 

? 
PORT SIDE TO APPROACH® ) 
STBD SIDE TO APPROACH ® ) 


S 
1 
S 
) 
> 
A 
A 


mid li COO Oh eos 
— « «>> VI<l UV) <i = 
ZOOM ewJtS SS eS Le 
MSYRUSZSZ dj ww SM LILI 
Ww mMoodatadtjwowodtd<ae 

COOtr  QOoOOkTRLMRMLOUOW 


E 
L 
T 


Pa 

m= ee OC 
= © © 
tC -— x 
LJ 
-— 


A 

9 

R 
PRPLO 


UL LL uJ LL 
cer we fa 
WilWWAcU = LULL 
LEewr Yet 
Wil WW LL 
aIqt <aerird 


Zt w 
<x<Mql «MEK X 
m— is DK ew —~ 
00 <m OI «KK O 
ZZmOZ wee 
Wiis 2 Lu << OO Ou 
QAA2VWAew 270. 
QAawaaduwa 
qIqoaaaqzwadc 

adr Wad 
SZz=—qt BQAagas 
Qogo =o2a CO 
CezOvtagsszaz 
Leck. Otte 

Ck See 
WWM MLOULiuwnm 
qq ni LL 
WartO@ZOW OFA 
=ODe = Le 
NoetxvrztIoOoz— 
YwCOHMNQIOwWM 

LJ =aAQ 
WWW O WLW MWI LL 7 
ode Praad 
pone Fd wt st ot 
me OF OOCOrFO 
2D eS DD ee es 
OO0OrF-OODFFFOZ 
CmUCLCOOUC- 
amazon 2ZzZZzawn 
cae Jems Sm [ite | aoe ome PP kaw pe Be 
NMLLOIWM LLL LL WI) 
Pd € 
<I fp = pf = J - 
COCKY eerrrxreo 
Qo em LU UU WU LU WL 
OGLNMNNMNMN nna 
ZRHeOSVS22Z222222z2Zz2™ 


LLE OZ) LLL pret pees etd fret fed Fed fend feng Peed 


S67 


50 


720 


509 


720 


52 


720 


50 


120 


INSERT TwO /7* CARDS HERE 





BIBLICG BAPHY 


Calvanc, C. N., An Investigation of the Stability of a 

System cf Iwo Ships Employing Automatic Control while 

On FParalleif Courses in Close fFroximity, 8.57 Thesis 

Hassachusetts Institute of Technology, cCambridce, 1970. 
rinciples of Naval Architecture, revised ed. 

a aeesuTe. The Society = Naval Architects and Marine 
ngineers, 1967. 

Pega tees of Naval Architecture, revised ed. 

ee t6, The Society of Naval Architects and Marine 
ngineéers, 1967. 

Naval Ship Research and Development Center Aol SEs 

202/67, Development of a Simulation Model of the. 

22S EBSEG AG SEG for a Destroyer Class Shig, by 
eL. Fatterson Jr., March 1969. 

Van, 1. V., Simulation of Propulsion Plant Dynawics and 

Their’ Effect cn Speed Control, 4.5. Thesis, Naval 

Postgraduate School, Monterey, 1974. 

Naval Ship Research and Development Center Report 

27-745, simulation of aaneuveriag Characteristics of a 

Desticyer Study Ship USing a Hodified Nenilinear Mcdel, 
y 5- H. Ercwn and k. Alvestad, August 197d. 


Thaler, G. J., "Computer Determination cf Low Order 
Mcdels fcr High Crder Systeas," Sn and 
Electrical Engineering, v. 2, p. TI7-123, T9775. 


66-5, The Analysis and Hodeling of Irregular Waves, by 
F. #. Sellars and T. A. Lcukakis, July 1966. 


Naval Shir Research and Development Center Report 








6-1C6y/7C, Digital Computer Simulation of Sea Wave 
Height ~ Time Series, by A. S. Fields, December S70. 
Lima, C. G., Multivariable Systems Design: A fwo Ship 

Controllier for Replenishment at Sea, H.S. Thesis, Naval 
ostcracuate School, Monterey, 1974. 

Astorquiza Vivar, G. M., Automatic Control System for 
Replenishment at Sea, M.S. Thesis, Naval Postgraduate 

Scacci, Hcnterey, 1975. 

Sarzetakis, T.,. Maneuvering Ccntrol of Replenishment at 
sey M.S. Thesis, Naval FEcstgraduate School, MNcnterey,~ 


Aguayo, E., Course Beepeng WieheAiicomtatic Control, M.S. 
Thesis, Naval Postgraduate School, Monterey, 1973. 


Mather, 1. S., Precision Surface Ship Separation 
tent and Maneuvering Control ee kad 
igués, paper presented at Ship Control Systé¢ 

Ui nd, Annapolis, Maryland, 4-5-6 Novemb 


Peers As Mancinwd, A. R., and Sarture, C. W., "A 


368 








16. 


17. 


18. 


19. 


20. 


26. 


Zl 


28. 


29. 


BIBLIOGRAPHY (cont.) 


Survey cf Adaptive Control Systems," I 
8 LMNs ere ontrol, v. PGAC-6, p. 10 


=m) Tran 


R nsactions 
2-108, December 


Leondes, C. f., Modern Control Systems Theory 
McGraw-Hill, $65. , 


Naticnal Aeronautics and Space Administration NASA 
CR-66S, Theoretical and Experimental Research on 
Digital Adaptive Control System, 5b eo Oe a Re Ge 
Barsh, &. &. Janitch, &. Re Chidag ara, and E. E. 
Buder, January 1967. 


Space Administration NASA 
heoretical and Experiment 
ive Control Systems, EY 
- &. Janitch and M. R. 


to 
tA otra 


Zaborszky, sh, 
Chidagrkara, July 19756. 

Naval Flectronics Laboratory NEL Report 1276, Adaptive 
Plant Identification Digital Control System, by W. V- 

Kershow, 30 March 196a. 


Defense Locumentation Center TR No. ASL-TR-61-27 Volume 
Fundéemental Study of Adaptive Control systems, bY Ree 


—_ a= a 


i 

E. Ralman, @. S. Englat, and BR. S. Bucy, April 
Naticnal Aeronautics and Space Administration NASA 
CR-715, Study of Opti gal and Adaptive Ccntrol fhecry, 
by C. Lb. Johnson, April 1967. . 


Mishkin, E. and Braun J 





Es 
Systems, McGraw-Hill, 1961. az —— 

Cadzcw, J. A., Martens, H. R., Discrete-Time and 
Comreter Control Systems, Prentice-Hall, Inc., 1973. 
Aguayc, E., Ccurse-Keeping with Automatic Control, M.S. 
T Ate val Bol coeadiace cha ie tcnser cram sise 
Hozos, A. Ge, Maneuvering Characteristics of 
Autcmatical Centroliied ships with Directio aity 
Unstable Hulis, 4.5. Thesis, Naval Postgraduate School, 
Bonterey, T974. 

Navy Marine Engineering Labcratory MEL Report 333/65 
Evaluation of Functional Performance of an Inteqrated 
Ship Contrcl Conning Console Dy Operator Personnel, by 
ne ae ne, W. J. Weingartner, and J. Cc. Townsend, 

ay : 


Hess, [L., "Shipboard saatend Reduction," All Hands, Nr. 
707, Fe 1421, December 1975. 


"Navy Ccmmands Report 4 More Ship Collisions," Navy 
Times, ¢~. 4, 31 December 1975. 


369 


INITIAL DISTRIBUTION LIST 


Defense Documentaticn Center 
Cameron Staticn 
Alexandria, Virginia 22314 


Litrary, Code 0212 
Naval Fcstgraduate School 
Monterey, Califcrnia 93940 


Cepartment Chairtan, Code 52 
Department of Electrical Engineering 
Naval Pcstgraduate School 

Mcnterey, California 93940 


Frofesscr George J. Thaler, Cede 52Tr 
Naval Pcstgraduate School 
Mcnterey, California 93940 


iPeeacineds Ubrin [IL, USN 
456 Linekerry Road 


Virginia Beach, Virginia 23452 


Réeidar Alvestad 
NSREIC Annapolis Lab 
Annagclis, Maryland 21402 


Asscc. Professor Alex Gerba, Jr., Code £2Gz 


Naval Pcstgraduate School 
Monterey, Califcrnia 93940 


Samuel E. Brown 
NSREC Annapolis Lab 
Annapclis, Maryland 21402 


370 





10. 


11. 


INITIAL DISTRIBUTION LIST (ccnt.) 


J. K. Zuidweg 1 
Royal Netherlands Naval College 
Den Helder, The Netherlands 


RADM BB. L. Walters 1 
Project Manager Surface Ship Froject Office PM18 


“Room 9SC8 National Center 3 


Department of the Navy 
Washirgtcn, DC 20362 
Attenticn PM18T 


Marvin CDCenicoff 1 
Office cf Naval Research Code 437 


MeNciaEgtcen,) Virginia 22217 








OPENS 





Sampled data adaptive 
diaital comouter contro) 
of surface ship maneuvers. 


Sampled data adaptive digital computer 


DUDLEY KNOX LIBRARY 





